BRANE- WORLD GRAVITY 



Roy Maartens 

Institute of Cosmology & Gravitation, Portsmouth University, Portsmouth POl 2EG, Britain 

The observable universe could be a 1+3-surface (the "brane") embedded in a 1+3+d-dimensional 
spacetime (the "bulk" ) , with Standard Model particles and fields trapped on the brane while gravity 
is free to access the bulk. At least one of the d extra spatial dimensions could be very large relative to 
the Planck scale, which lowers the fundamental gravity scale, possibly even down to the electroweak 
(~TeV) level. This revolutionary picture arises in the framework of recent developments in M theory. 
The 1+10-dimensional M theory encompasses the known 1+9-dimensional superstring theories, and 
is widely considered to be a promising potential route to quantum gravity. General relativity cannot 
describe gravity at high enough energies and must be replaced by a quantum gravity theory, picking 
up significant corrections as the fundamental energy scale is approached. At low energies, gravity 
is localized at the brane and general relativity is recovered, but at high energies gravity "leaks" 
into the bulk, behaving in a truly higher-dimensional way. This introduces significant changes to 
gravitational dynamics and perturbations, with interesting and potentially testable implications for 
high-energy astrophysics, black holes and cosmology. Brane-world models offer a phenomenological 
way to test some of the novel predictions and corrections to general relativity that are implied by 
M theory. This review discusses the geometry, dynamics and perturbations of simple brane-world 
models for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane- worlds based 
on the Randall-Sundrum models. 



At high enough energies, Einstein's theory of general relativity breaks down, and will be superceded by a quantum 
gravity theory. The classical singularities predicted by general relativity in gravitational collapse and in the hot big 
bang will be removed by quantum gravity. But even below the fundamental energy scale that marks the transition to 
quantum gravity, significant corrections to general relativity will arise. These corrections could have a major impact 
on the behaviour of gravitational collapse, black holes and the early universe, and they could leave a trace-a "smoking 
gun" -in various observations and experiments. Thus it is important to estimate these corrections and develop tests for 
detecting them or ruling them out. In this way, quantum gravity can begin to be subject to testing by astrophysical 
and cosmological observations. 

Developing a quantum theory of gravity and a unified theory of all the forces and particles of nature are the two main 
goals of current work in fundamental physics. There is as yet no generally accepted (pre-)quantum gravity theory. 
Two of the main contenders are M theory and quantum geometry (loop quantum gravity) 0- It is important to 
explore the astrophysical and cosmological predictions of both these approaches. This review considers only models 
that arise within the framework of M theory, and mainly the 5-dimensional warped brane- worlds. 



One of the fundamental aspects of string theory is the need for extra spatial dimensions. This revives the original 
higher-dimensional ideas of Kaluza and Klein in the 1920's, but in a new context of quantum gravity. An important 
consequence of extra dimensions is that the 4-dimensional Planck scale M p = M4 is no longer the fundamental scale, 
which is Mi + d, where d is the number of extra dimensions. This can be seen from the modification of the gravitational 
potential. For an Einstein-Hilbert gravitational action we have, 



I. INTRODUCTION 



A. Heuristics of higher-dimensional gravity 
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The static weak field limit of the field equations leads to the 4 + rf-dimensional Poisson equation, whose solution is 
the gravitational potential, 

V(r) oc jfg . (1.4) 

If the length scale of the extra dimensions is L, then on scales r < L, the potential is 4 + rf-dimensional, V ~ r~( 1+d \ 
By contrast, on scales large relative to L, where the extra dimensions do not contribute to variations in the potential, 
V behaves like a 4-dimensional potential, i.e., r ~ L in the d extra dimensions, and V ~ L~ d r~ 1 . This means that 
the usual Planck scale becomes an effective coupling constant, describing gravity on scales much larger than the extra 
dimensions, and related to the fundamental scale via the volume of the extra dimensions: 

Ml~MH d d L d . (1.5) 

If the extra-dimensional volume is Planck scale, i.e. L ~ M~ , then Mj+d ~ M p . But if the extra-dimensional volume 
is significantly above Planck-scale, then the true fundamental scale Mi +d can be much less than the effective scale 
M p ~ 10 19 GeV. In this case, we understand the weakness of gravity as due to the fact that it "spreads" into extra 
dimensions and only a part of it is felt in 4 dimensions. 

A lower limit on Mi +d is given by null results in table-top experiments to test for deviations from Newton's law in 
4 dimensions, V oc r~ 1 . These experiments currently Q probe sub-millimetre scales, so that 

L < KT 1 mm ~ (1(T 15 TeV)" 1 M 4+d > i ( 32 - 15d )/( d + 2 ) TeV . (1.6) 

Stronger bounds for brane-worlds with compact flat extra dimensions can be derived from null results in particle 
accelerators and in high-energy astrophysics 
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B. Brane-worlds and M theory 



String theory thus incorporates the possibility that the fundamental scale is much less than the Planck scale felt 
in 4 dimensions. There are five distinct 1+9-dimensional superstring theories, all giving quantum theories of gravity. 
Discoveries in the mid-90s of duality transformations that relate these superstring theories and the 1+10-dimensional 
supergravity theory, led to the conjecture that all of these theories arise as different limits of a single theory, which 
has come to be known as M theory. The 11th dimension in M theory is related to the string coupling strength; the 
size of this dimension grows as the coupling becomes strong. At low energies, M theory can be approximated by 
1+10-dimensional supergravity. 

It was also discovered that p-branes, which are extended objects of higher dimension than strings (1-branes), play 
a fundamental role in the theory. In the weak coupling limit, p-branes (p > 1) become infinitely heavy, so that they 
do not appear in the perturbative theory. Of particular importance among p-branes are the D-branes, on which open 
strings can end. Roughly speaking, open strings, which describe the non-gravitational sector, are attached at their 
endpoints to branes, while the closed strings of the gravitational sector can move freely in the bulk. Classically, this 
is realised via the localization of matter and radiation fields on the brane, with gravity propagating in the bulk (see 
Fig. 1). 

In the Horava-Witten solution pj, gauge fields of the standard model are confined on two 1+9-branes located at 
the end points of an S 1 /Z% orbifold, i.e., a circle folded on itself across a diameter. The 6 extra dimensions on the 
branes are compactified on a very small scale, close to the fundamental scale, and their effect on the dynamics is felt 
through "moduli" fields, i.e. 5D scalar fields. A 5D realization of the Horava-Witten theory and the corresponding 
brane- world cosmology is given in 

These solutions can be thought of as effectively 5-dimensional, with an extra dimension that can be large relative to 
the fundamental scale. They provide the basis for the Randall-Sundrum 2-brane models of 5-dimensional gravity 
(see Fig. 2). The single-brane Randall-Sundrum models 10] with infinite extra dimension arise when the orbifold 
radius tends to infinity. The RS models are not the only phenomenological realizations of M theory ideas. They were 
preceded by the Arkani-Hamed-Dimopoulos-Dvali (ADD) brane-world models, which put forward the idea that 
a large volume for the compact extra dimensions would lower the fundamental Planck scale, 

M ew ~ 1 TeV < M A+d < M p ~ 10 16 TeV , (1.7) 

where M GW is the electroweak scale. If M^ +d is close to the lower limit in Eq. 11.71) . then this would address the 
long-standing "hierarchy" problem, i.e. why there is such a large gap between M ew and M p . 
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FIG. 1: Schematic of confinement of matter to the brane, while gravity propagates in the bulk (from 



In the ADD models, more than one extra dimension is required for agreement with experiments, and there is 
"democracy" amongst the equivalent extra dimensions, which are typically flat. By contrast, the RS models have 
a "preferred" extra dimension, with other extra dimensions treated as ignorable (i.e., stabilized except at energies 
near the fundamental scale). Furthermore, this extra dimension is curved or "warped" rather than flat: the bulk is a 
portion of anti de Sitter (AdSs) spacetime. As in the Horava-Witten solutions, the RS branes are ^-symmetric (mirror 
symmetry), and have a tension, which serves to counter the influence of the negative bulk cosmological constant on 
the brane. This also means that the self-gravity of the branes is incorporated in the RS models. The novel feature 
of the RS models compared to previous higher-dimensional models is that the observable 3 dimensions are protected 
from the large extra dimension (at low energies) by curvature rather than straightforward compactification. 

The RS brane- worlds and their generalizations (to include matter on the brane, scalar fields in the bulk, etc.) provide 
phenomenological models that reflect at least some of the features of M theory, and that bring exciting new geometric 
and particle physics ideas into play. The RS2 models also provide a framework for exploring holographic ideas that 
have emerged in M theory. Roughly speaking, holography suggests that higher-dimensional gravitational dynamics 
may be determined from knowledge of the fields on a lower-dimensional boundary. The AdS/CFT correspondence 
is an example, in which the classical dynamics of the higher-dimensional gravitational field are equivalent to the 
quantum dynamics of a conformal field theory (CFT) on the boundary. The RS2 model with its AdSs metric satisfies 
this correspondence to lowest perturbative order ^| ( see a ls° for the AdS / CFT correspondence in a cosmological 
context) . 

In this review, I focus on RS brane-worlds (mainly RS 1-brane) and their generalizations, with the emphasis on 
geometry and gravitational dynamics (see |14l lla. Il6l for previous reviews with a broadly similar approach). Other 
recent reviews focus on string-theory aspects, e.g. [jj]], or on particle physics aspects, e.g. Before turning to a 

more detailed analysis of RS brane-worlds, I discuss the notion of Kaluza-Klein (KK) modes of the graviton. 



The dilution of gravity via extra dimensions not only weakens gravity on the brane, it also extends the range of 
graviton modes felt on the brane beyond the massless mode of 4-dimensional gravity. For simplicity, consider a flat 
brane with one flat extra dimension, compactified through the identification y «-> y + 2irnL, where n = 0, 1, 2, • • • . 
The perturbative 5D graviton amplitude can be Fourier expanded as 



C. 



Heuristics of KK modes 




(1.8) 



n 



where /„ are the amplitudes of the KK modes, i.e. the effective 4D modes of the the 5D graviton. To see that these 
KK modes are massive from the brane viewpoint, we start from the 5D wave equation that the massless 5D field / 
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FIG. 2: The RS 2-brane model (from @). 



satisfies (in a suitable gauge): 

( 5 ta/ = o □/ + a y 2 / = o. (i.9) 

It follows that the KK modes satisfy a 4D Klein-Gordon equation with an effective 4D mass, m n , 

n/n = ^ /n , m n = — . (1.10) 

The massless mode, /o, is the usual 4D graviton mode. But there is a tower of massive modes, L ,2L , ■ ■ ■ , 
which imprint the effect of the 5D gravitational field on the 4D brane. Compactness of the extra dimension leads to 
discreteness of the spectrum. For an infinite extra dimension, L — > oo, the separation between the modes disappears 
and the tower forms a continuous spectrum. In this case, the coupling of the KK modes to matter must be very weak 
in order to avoid exciting the lightest massive modes with m > 0. 

From a geometric viewpoint, the KK modes can also be understood via the fact that the projection of the null 
graviton 5-momentum ^pa onto the brane is timelike. If the unit normal to the brane is ua, then the induced metric 
on the brane is 

9AB = {5) 9AB ~ n A n B , {5) g A Bn A n B = 1 , g AB n B = , (1.11) 

and the 5-momentum may be decomposed as 

{ % A = mriA+PA , PAn A = 0, m = (5) p A n A , (1.12) 

where pa — gAB ^p B is the projection along the brane, depending on the orientation of the 5-momentum relative to 
the brane. The effective 4-momentum of the 5D graviton is thus pa- Expanding = 0, we find that 

g A BP A p B = -m 2 - (1.13) 

It follows that the 5D graviton has an effective mass m on the brane. The usual 4D graviton corresponds to the zero 
mode, m — 0, when ^pa is tangent to the brane. 

The extra dimensions lead to new scalar and vector degrees of freedom on the brane. In 5D, the spin-2 graviton is 
represented by a metric perturbation ^Iiab that is transverse traceless: 

^h A A = = d B ^h A B ■ (1-14) 
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In a suitable gauge, ^\ab contains a 3D transverse traceless perturbation hij, a 3D transverse vector perturbation 
Sj and a scalar perturbation (3, which each satisfy the 5D wave equation |19|: 

( ^h AB — » hij, (1.15) 

= o = a,/^' , = 0, (i.i6) 

a 

(D + d 2 ) I |=0. (1.17) 

hi. 

The other components of ^Hab are determined via constraints once these wave equations are solved. The 5 degrees 
of freedom (polarizations) in the 5D graviton are thus split into 2 (hij) + 2 (£j) +1 (fi) degrees of freedom in 4D. On 
the brane, the 5D graviton field is felt as 

• a 4D spin-2 graviton hij (2 polarizations) 

• a 4D spin-1 gravi-vector (gravi-photon) Sj (2 polarizations) 

• a 4D spin-0 gravi-scalar j3. 

The massive modes of the 5D graviton are represented via massive modes in all 3 of these fields on the brane. The 
standard 4D graviton corresponds to the massless zero- mode of hij. 

In the general case of d extra dimensions, the number of degrees of freedom in the graviton follows from the 
irreducible tensor representations of the isometry group as ^ (d + 1) (d + 4) . 

II. RANDALL-SUNDRUM BRANE- WORLDS 

RS brane-worlds do not rely on compactification to localize gravity at the brane, but on the curvature of the bulk 
(sometimes called "warped compactification"). What prevents gravity from 'leaking' into the extra dimension at low 
energies is a negative bulk cosmological constant, 

A 5 = -| = -6/i 2 , (2.1) 

where £ is the curvature radius of AdSs and /i is the corresponding energy scale. The curvature radius determines the 
magnitude of the Riemann tensor: 



{5) Rabcd - —p f 



(5 W {5) g B D - {5) gAD (5 W 



(2.2) 



The bulk cosmological constant acts to "squeeze" the gravitational field closer to the brane. We can see this clearly 
in Gaussian normal coordinates X A = (a;' 1 , y) based on the brane at y = 0, for which the AdSs metric takes the form 

( 5 fe 2 = e- 2 ^'% v dx tx dx v + dy 2 , (2.3) 

with rjfn, the Minkowski metric. The exponential warp factor reflects the confining role of the bulk cosmological 
constant. The ^-symmetry about the brane at y — is incorporated via the \y\ term. In the bulk, this metric is a 
solution of the 5D Einstein equations, 

(5) Gab = -A 5 (5 W , (2-4) 

i.e., ^Tab = in Eq. (|1.2fl . The brane is a flat Minkowski spacetime, gAB(x a ,0) = rj^S 11 a8 u b, with self-gravity in 
the form of brane tension. One can also use Poincare coordinates, which bring the metric into manifestly conformally 
flat form: 

(2 

( 5 )ds 2 = — [ij^dx^dx" + dz 2 ] , (2.5) 



where z = £e v ^ e . 

The two RS models are distinguished as follows: 
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RS 2-BRANE: 

There are two branes in this model at y = and y = L, with .^-symmetry identifications 

V + L^L-y. (2.6) 
The branes have equal and opposite tensions, ±A, where 

The positive-tension brane has fundamental scale M5 and is "hidden". Standard model fields are confined on 
the negative tension (or "visible" ) brane. Because of the exponential warping factor, the effective scale on the 
visible brane at y = L is M p , where 



Mp = Ml I 



1 - e 



-2L/1 



(2.8) 



So the RS 2-brane model gives a new approach to the hierarchy problem. Because of the finite separation 
between the branes, the KK spectrum is discrete. Furthermore, at low energies gravity on the branes becomes 
Brans-Dicke-like, with the sign of the Brans-Dicke parameter equal to the sign of the brane tension 20] . In order 
to recover 4D general relativity at low energies, a mechanism is required to stabilize the inter-brane distance, 
which corresponds to a scalar field degree of freedom known as the radion 0, [2^ . 

RS 1-BRANE: 

In this model [lOj, there is only one, positive tension, brane. It may be thought of as arising from sending the 
negative tension brane off to infinity, L — > 00. Then the energy scales are related via 

a Ml 

Ml = -f . (2.9) 
The infinite extra dimension makes a finite contribution to the 5D volume because of the warp factor: 

J d 5 X^Wg = 2 J d 4 x J™ dye-^l l = l - J d 4 x. (2.10) 
Thus the effective size of the extra dimension probed by the 5D graviton is I. 



I will concentrate mainly on RS 1-brane from now on, referring to RS 2-brane occasionally. The RS 1-brane 
models are in some sense the most simple and geometrically appealing form of brane-world model, while at the same 
time providing a framework for the AdS/CFT correspondence [l2|,[RJ- The RS 2-brane models introduce the added 
complication of radion stabilization, as well as possible complications arising from negative tension. However, they 
remain important and will occasionally be discussed. 

In RS 1-brane, the negative A5 is offset by the positive brane tension A. The fine-tuning in Eq. (|2.7|l ensures 
that there is zero effective cosmological constant on the brane, so that the brane has the induced geometry of 
Minkowski spacetime. To see how gravity is localized at low energies, we consider the 5D graviton perturbations of 
the metric [hJ HJ E|: 

{5} 9ab - (5 W + e- 2 ^ e ^h AB , ^h Ay = = = (%% . (2.11) 

(See Fig. 3.) This is the RS gauge, which is different from the gauge used in Eq. i|1.15fl . but which also has no 
remaining gauge freedom. The 5 polarizations of the 5D graviton are contained in the 5 independent components of 
^h^ u in the RS gauge. 

We split the amplitude / of ^ 5 'Hab m to 3D Fourier modes, and the linearized 5D Einstein equations lead to the 
wave equation (y > 0) 



f + k 2 f = e 



-2y/i 



r-\r 



(2.12) 
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FIG. 3: Gravitational field of small point particle on the brane in RS gauge (from 20]). 



Separability means we can write 



and the wave equation reduces to 



(p m + (m 2 + k 2 )(p rn = , 
4 
1 



The zero mode solution is 



<po(t) = A 0+ e +ikt + A _e- m , 
fo(y) = B Q + C e 4 y/ e , 

and the m > solutions are 

ip m (t) = A m+ exp m 2 + k 2 tj + A m - cxp i\j m 2 + k 2 tj , 

f m (y) = e 2y / e [B m J 2 (mle^ + C m Y 2 {mie v l 1 } 



(2.13) 

(2.14) 
(2.15) 



(2.16) 
(2.17) 



(2.18) 
(2.19) 



The boundary condition for the perturbations arises from the junction conditions, Eq. 13.19|) . discussed below, 
and leads to f'(t,0) = 0, since the transverse traceless part of the perturbed energy- momentum tensor on the brane 
vanishes. This implies 



Cq — , C„ 



Jx{ml) 
Yi(m£) 



B n 



The zero mode is normalizable, since 



/>OG 

/ B Q e-^% 
Jo 



< oo . 



(2.20) 



(2.21) 
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Its contribution to the gravitational potential V — i ' 5 Vioo gives the 4D result, V oc r 1 . The contribution of the 
massive KK modes sums to a correction of the 4D potential. For r <C £, one obtains 

V(r)-^«. (2-22, 
which simply reflects the fact that the potential becomes truly 5D on small scales. For r £, 



GM A . 2£ 2 
3r 



V(r) « — ( 1 + — ) , (2.23) 



which gives the small correction to 4D gravity at low energies from extra-dimensional effects. These effects serve to 
slightly strengthen the gravitational field, as expected. 

Table-top tests of Newton's laws currently find no deviations down to (^(lCP 1 mm), so that £ < 0.1 mm in Eq. I|2.23[l . 
Then by Eqs. (|2.7(l and l|2.9|l . this leads to lower limits on the brane tension and the fundamental scale of the RS 
1-brane model: 

A > (1 TeV) 4 , M 5 > 10 5 TeV. (2.24) 

These limits do not apply to the 2-brane case. 

For the 1-brane model, the boundary condition, Eq. Q2.20[l . admits a continuous spectrum m > of KK modes. In 
the 2-brane model, f'(t,L) — must hold in addition to Eq. (|2.2l)[l . This leads to conditions on m, so that the KK 
spectrum is discrete: 

m„ = y e~ L / £ where Y x (m£)Ji{x n ) = Jx(m£)Y x {x n ) . (2.25) 

The limit Eq. (|2.24|) indicates that there are no observable collider, i.e. O(TeV), signatures for the RS 1-brane model. 
The 2-brane model by contrast, for suitable choice of L and £ so that m x — O(TeV), does predict collider signatures 
that are distinct from those of the ADD models . 



III. COVARIANT APPROACH TO BRANE- WORLD GEOMETRY AND DYNAMICS 

The RS models and the subsequent generalization from a Minkowski brane to a Friedmann-Robertson- Walker 
(FRW) brane 131111113; were derived as solutions in particular coordinates of the 5D Einstein equations, together 
with the junction conditions at the ^-symmetric brane. A broader perspective, with useful insights into the inter- 
play between 4D and 5D effects, can be obtained via the covariant Shiromizu-Maeda-Sasaki approach 27], in which 
the brane and bulk metrics remain general. The basic idea is to use the Gauss-Codazzi equations to project the 5D 
curvature along the brane. (The gen eral formalism for relating the geometries of a spacetime and of hypersurfaces 
within that spacetime is given in |28|.) 

The 5D field equations determine the 5D curvature tensor; in the bulk, they are 

^G AB = -A s ( 5 W + 4 (5) Tab , (3.1) 

where ^Tab represents any 5D energy-momentum of the gravitational sector (e.g., dilaton and moduli scalar fields, 
form fields). 

Let y be a Gaussian normal coordinate orthogonal to the brane (which is at y = without loss of generality), so 
that nAdX A = dy, with n A the unit normal. The 5D metric in terms of the induced metric on {y = const} surfaces 
is locally given by 

(5 W = 9AB + n A n B , ^ds 2 = g^{x a ,y)dx tl dx v + dy 2 . (3.2) 

The extrinsic curvature of {y = const} surfaces describes the embedding of these surfaces. It can be defined via the 
Lie derivative or via the covariant derivative: 

K AB = \ £„ 9ab - 9A C (5 V c «i3 , (3.3) 

so that 



K [AB] = = K AB n B , (3.4) 
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where square brackets denote anti-symmetrization. The Gauss equation gives the 4D curvature tensor in terms of the 
projection of the 5D curvature, with extrinsic curvature corrections: 

Rabcd = ^'Refgh 9a E 9b F gc G 9d H + 2K A [ C K D ] B , (3.5) 

and the Codazzi equation determines the change of K AB along {y = const} via 

V B K B A ~ W A K = ^R BC g A B n c , (3.6) 

where K = K A A . 

Some other useful projections of the 5D curvature are: 

(5) Refgh 9A E 9B F 9c G n H = 2V [A K B]C , (3.7) 

(5) Refgh 9A E n F g B G n H = -£ n K AB + K AC K C B , (3.8) 

(5) Rcd 9a C 9b D = Rab - £ n K AB - K K AB + 2K AC K C B . (3.9) 

The 5D curvature tensor has Weyl (trace-free) and Ricci parts: 

{5) Rabcd = (5) C AC bd + \ { (5 W[c (5) i?D]B - i5) 9B[c (5) R D] a} - W* W (5)i? ■ (3.10) 



A. Field equations on the brane 

Using Eqs. (|3.1|l and (|3.5|) . it follows that 

1 . 2 



Gnu 



3 5 



{5} TABg» A g„ B + {5) T AB n A n B 



,(5>T 



9nv 



KK^ - K^K ay + i [K^K a(i - K 2 ] g^ - £, 



where ^ 5 ^T = ^ 5 'T A A , and where 

£ MV = {5) C AC BDn c n D g^ A gv B , 
is the projection of the bulk Weyl tensor orthogonal to n A . This tensor satisfies 

£abit- B = = £[ab] = £a A , 



(3.11) 



(3.12) 



(3.13) 



by virtue of the Weyl tensor symmetries. Evaluating Eq. I|3.11|l on the brane (strictly, as y — > ±0, since £ AB is not 
defined on the brane |2]j) will give the field equations on the brane. 

First, we need to determine K^ v at the brane from the junction conditions. The total energy-momentum tensor on 
the brane is 



rpb 



= Tuu - Ao 



(3.14) 



where is the energy-momentum tensor of particles and fields confined to the brane (so that T AB n B = 0). The 5D 
field equations, including explicitly the contribution of the brane, are then 



(5) G 



AB 



-A,( 5 ), 



5 v yAB + K 5 



(5>r 



AB 



T 



brane 
AB 



S(y) 



(3.15) 



Here the delta function enforces in the classical theory the string theory idea that Standard Model fields are confined 
to the brane. This is not a gravitational confinement, since there is in general a nonzero acceleration of particles 
normal to the brane [35| . 

Integrating Eq. (|3.15|) along the extra dimension from y = — e to y = +e, and taking the limit e — * 0, leads to the 
Israel-Darmois junction conditions at the brane, 



9 tut g 



o. 



Kt„ - K„„ = -k 



T, 



3 



(3.16) 
(3.17) 
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where T brane = g fiU T^ ane . The Z2 symmetry means that when you approach the brane from one side and go through 
it, you emerge into a bulk that looks the same, but with the normal reversed, n A — > —n A . Then Eq. I|3.3|) implies 
that 



K~ = -K + 



(3.18) 



so that we can use the junction condition Eq. I|3.17|l to determine the extrinsic curvature on the brane: 



K 



~ 2 



T ^ + - (A - T) 



(3.19) 



where T = T%, we have dropped the (+) and we evaluate quantities on the brane by taking the limit y — > +0. 
Finally we arrive at the induced field equations on the brane, by substituting Eq. I|3.19[) into Eq. I|3.11|l : 



2 K 2 K 2 



(3.20) 



The 4D gravitational constant is an effective coupling constant inherited from the fundamental coupling constant, 
and the 4D cosmological constant is nonzero when the RS balance between the bulk cosmological constant and the 
brane tension is broken: 



2 _ 2 _ 1 , 

r€ — f^A — „ y\Ki 



5 > 



A = - [A 5 + k 2 X] , 



(3.21) 
(3.22) 



The first correction term relative to Einstein's theory is quadratic in the energy-momentum tensor, arising from 
the extrinsic curvature terms in the projected Einstein tensor: 



Spa* — ~[2 TTt * u ~ ^T^aT",, + [STapT 0113 - T 2 ] 



(3.23) 



The second correction term is the projected Weyl term. The last correction term on the right of Eq. (|3.20() . which 
generalizes the field equations in [27|. is 



A B . 
AB9p, 9v + 



^>T AB n A n B - - 



9y.v 



(3.24) 



where ^ 5 ^T AB describes any stresses in the bulk apart from the cosmological constant (see |2j| for the case of a scalar 
field). 

What about the conservation equations? Using Eqs. (|3.1() . I|3.6|) and l|3.19[) . one obtains 



V%„ = -2 {5) T AB n A g B ^ 



(3.25) 



Thus in general there is exchange of energy-momentum between the bulk and the brane. From now on, we will assume 
that 



AB 



=> T^ = 0. 



so that 



^G AB = -A 5 <®g AB , (in the bulk) 



M „ - -A.9, w + K^Tpu + 6— S M „ ~ E^ v , (on the brane) 



and one then recovers from Eq. 13.25[) the standard 4D conservation equations, 

vt^ = . 

This means that there is no exchange of energy-momentum between the bulk and the brane; their interaction is purely 
gravitational. Then the 4D contracted Bianchi identities (VG^ = 0), applied to Eq. I|3.20|l . lead to 



(3.26) 

(3.27) 
(3.28) 

(3.29) 



6k 2 



(3.30) 



which shows qualitatively how 1+3 spacetime variations in the matter-radiation on the brane can source KK modes. 

The induced field equations l|3.28|) show two key modifications to the standard 4D Einstein field equations arising 
from extra-dimensional effects: 
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• ~ {Tfiv} 2 is the high-energy correction term, which is negligible for p <C A, but dominant for p ^> A: 



|k 2 VI a a 



(3.31) 



• Ef+v, the projection of the bulk Weyl tensor on the brane, encodes corrections from 5D graviton effects (the KK 
modes in the linearized case). From the brane-observer viewpoint, the energy-momentum corrections in S^ v are 
local, whereas the KK corrections in E ^ v are nonlocal, since they incorporate 5D gravity wave modes. These 
nonlocal corrections cannot be determined purely from data on the brane. In the perturbative analysis of RS 
1-brane which leads to the corrections in the gravitational potential, Eq. (|2.23|) . the KK modes that generate 
this correction are responsible for a nonzero S^ v \ this term is what carries the modification to the weak-field 
field equations. The 9 independent components in the tracefree £ M „ are reduced to 5 degrees of freedom by 
Eq. I|3.30J1 ; these arise from the 5 polarizations of the 5D graviton. 

Note that the covariant formalism applies also to the two-brane case. In that case, the gravitational influence 
of the second brane is felt via its contribution to £^ u . 



B. 5-dimensional equations and the initial-value problem 



The effective field equations are not a closed system. One needs to supplement them by 5D equatio ns g overning 



£fj,i,, which are obtained from the 5D Einstein and Bianchi equations. This leads to the coupled system |3 

1 



G 



1 



£ n ^w = V Q S Q ( Ml/ ) + -A 5 {K^ — g^ v K) + K a ^R fi aup 



(fj,£i/)a 



KE 



F E 



+ (K„ a K v p - K aP K^) K al3 , 

where the "magnetic" part of the bulk Weyl tensor, counterpart to the "electric" part £ M „, is 

Bfj, ua = 9n A 9v B g a C {5) C ABC Dn D . 
These equations are to be solved subject to the boundary conditions at the brane, 



B,j, va = 2V[^K v ] a = — «iV[^T„] a - -gu]aTj; 



(3.32) 



(3.33) 
(3.34) 
(3.35) 



(3.36) 

(3.37) 
(3.38) 



where A = B denotes A | brane = -B| brane- 

The above equations have been used to develop a covariant analysis of the weak field [30( . They can also be used 
to develop a Taylor expansion of the metric about the brane. In Gaussian normal coordinates, Eq. (|3.2|l . we have 
£n = d/dy. Then we find 



g^u{x,y) = g^ v (x,0) - k 5 
1 



—£fMu + 7^5 ^ T fla T a u + -(A - T)T Miy 



y=o+ 

-{- K U\-Tf 
6 I 6 5V ; 



A, 



(3.39) 



y=o+ 



In a non-covariant approach based on a specific form of the bulk metric in particular coordinates, the 5D Bianchi 
identities would be avoided and the equivalent problem would be one of solving the 5D field equations, subject to 
suitable 5D initial conditions and to the boundary conditions Eq. (|3.19|) on the metric. The advantage of the covariant 
splitting of the field equations and Bianchi identities along and normal to the brane is the clear insight that it gives 
into the interplay between the 4D and 5D gravitational fields. The disadvantage is that the splitting is not well suited 
to dynamical evolution of the equations. Evolution off the timelike brane in the spacelike normal direction does not 
in general constitute a well-defined initial value problem |3lj . One needs to specify initial data on a 4D spacelike (or 
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null) surface, with boundary conditions at the brane(s) ensuring a consistent evolution [33]. Clearly the evolution of 
the observed universe is dependent upon initial conditions which are inaccessible to brane-bound observers; this is 
simply another aspect of the fact that the brane dynamics is not determined by 4D but by 5D equations. The initial 
conditions on a 4D surface could arise from models for creation of the 5D universe ptj, l34| . from dynamical attractor 
behaviour 33] or from suitable conditions (such as no incoming gravitational radiation) at the past Cauchy horizon 
if the bulk is asymptotically AdS. 



C. The brane viewpoint: a 1+3-covariant analysis 

A systematic analysis can be developed from the viewpoint of a brane-bound observer, following j^. The effects 
of bulk gravity are conveyed, from a brane observer viewpoint, via the local (S^) and nonlocal (S^) corrections to 
Einstein's equations. (In the more general case, bulk effects on the brane are also carried by which describes 
any 5D fields.) The term cannot in general be determined from data on the brane, and the 5D equations above 
(or their equivalent) need to be solved in order to find £a V . 

The general form of the brane energy-momentum tensor for any matter fields (scalar fields, perfect fluids, kinetic 
gases, dissipative fluids, etc.), including a combination of different fields, can be covariantly given in terms of a chosen 
4-velocity it M as 

T M „ = pu^u u + ph^ + n^v + q^u v + q^u^ . (3.40) 
Here p and p are the energy density and isotropic pressure, and 

hpv = g^u + UftUv = (5 W<^ - n v n v + u n u » > ( 3 - 41 ) 

projects into the comoving rest space orthogonal to it M on the brane. The momentum density and anisotropic stress 
obey 

9m = <?<m} ' = ^(nv) > (3.42) 
where angled brackets denote the spatially projected, symmetric and tracefree part: 



W a0 . (3.43) 



v w = h lt v v„, w (liv) = 

In an inertial frame at any point on the brane, we have 

u" = (1, 0) , V = diag(0, 1, 1, 1) , V» = (0, V t ) , W M o = = ^ W u = - W jt , (3.44) 
where i, j — 1, 2, 3. 

The tensor S^, which carries local bulk effects onto the brane, may then be irreducibly decomposed as 

- j^(p + 2p)7r M „ + ■n a (^y) a + q^u) + ^pq(^u u ) - -^q a n a ^u u ) . (3.45) 
This simplifies for a perfect fluid or minimally-coupled scalar field: 



Sfj,u = j^P \pu^u u + (p + 2p) h^] . (3.46) 



The trace free carries nonlocal bulk effects onto the brane, and contributes an effective "dark" radiative energy- 
momentum on the brane, with energy density pz, pressure pg/3, momentum density q^ and anisotropic stress ir^: 



-^Z^u = Ps (u^Uu + ) + q^u v + q^Ufi + . (3.47) 



We can think of this as a KK or Weyl "fluid" . The brane "feels" the bulk gravitational field through this effective 
fluid. More specifically: 
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• The KK (or Weyl) anisotropic stress ir^ v incorporates the scalar or spin-0 ("Coulomb"), the vector (transverse) 
or spin-1 (gravimagnetic) and the tensor (transverse traceless) or spin-2 (gravitational wave) 4D modes of the 
spin-2 5D graviton. 

• The KK momentum density incorporates spin-0 and spin-1 modes, and defines a velocity of the Weyl fluid 
relative to m m via = pgv^. 

• The KK energy density ps , often called the "dark radiation" , incorporates the spin-0 mode. 

In special cases, symmetry will impose simplifications on this tensor. For example, it must vanish for a conformally 
flat bulk, including AdSs, 

(5 Wb conformally flat =>■ = . (3.48) 

The RS models have a Minkowski brane in an AdSs bulk. This bulk is also compatible with an FRW brane. However, 
the most general vacuum bulk with a Friedmann brane is Schwarzschild-anti de Sitter spacetime |36| . Then it follows 
from the FRW symmetries that 

Schwarzschild AdS 5 bulk, FRW brane: cfc=a = rf v , (3.49) 

where ps = only if the mass of the black hole in the bulk is zero. The presence of the bulk black hole generates via 
Coulomb effects the dark radiation on the brane. 

For a static spherically symmetric brane (e.g. the exterior of a static star or black hole) |37) . 

static spherical brane: q^ = . (3.50) 

This condition also holds for a Bianchi I brane j3^. In these cases, 7rf„ is not determined by the symmetries, but by 
the 5D field equations. By contrast, the symmetries of a Godel brane fix |39| . 

The brane-world corrections can conveniently be consolidated into an effective total energy density, pressure, mo- 
mentum density and anisotropic stress. 

Ptot = p + (2p 2 - tor^if) + p £ , (3.51) 

Ptot - P + (2p 2 + 4pp + tt^tt^ - Aq^) + y , (3.52) 

C = «/* + ^ ( 2 P9/* - + 4 , (3-53) 



7T„„ = 7T 



nv - »»v t i~(p + 3p)w^ v + 3w a ( li ir 1/ ) a + 3g <At g v )J + . (3.54) 

These general expressions simplify in the case of a perfect fluid (or minimally coupled scalar field, or isotropic one- 
particle distribution function), i.e., for q^ = = 7ty„: 

Ptot - P+ {x {2P + P) + ! J ' (3 ' 56) 

C = 4 , (3-57) 

^ = rfv ■ (3-58) 

Note that nonlocal bulk effects can contribute to effective imperfect fluid terms even when the matter on the brane 
has perfect fluid form: there is in general an effective momentum density and anisotropic stress induced on the brane 
by massive KK modes of the 5D graviton. 

The effective total equation of state and sound speed follow from Eqs. (|3.55|) and l|3.5(j|l as 

p tot w + {l + 2w)p/2\ + p £ /3p 

Wtot = = 7- 777T j ] , 3.59) 

Ptot l + p 2\ + p £ p 



„2 _ Ptot 
-tot — • 

Ptot 



2 , P + P 4 P£ 



p+X 9(p + p)(l + p/X) 



4p £ 



3(p + p)(l + p/A) 



(3.60) 
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where w = p/p and cj? = p/p. At very high energies, i.e., p > A, we can generally neglect ps (e.g., in an inflating 
cosmology), and the effective equation of state and sound speed are stiffened: 

iwtot ~ 2w + 1 , Ct ot ps + w + 1 . (3.61) 

This can have important consequences in the early universe and during gravitational collapse. For example, in a very 
high-energy radiation era, w = 5, the effective cosmological equation of state is ultra-stiff: Wtot ~ f ■ I n late-stage 
gravitational collapse of pressureless matter, w = 0, the effective equation of state is stiff, u> to t ~ 1, and the effective 
pressure is nonzero and dynamically important. 



D. Conservation equations 



Conservation of gives the standard general relativity energy and momentum conservation equations, in the 
general, nonlinear case: 

p + 6(p + p) + V% + 2A% + a^Tr^ = , (3.62) 

4 

4<m> + 3@9p + V M p + (p + p)A M + + + a^q" - e^ a u»q a = . (3.63) 

In these equations, an overdot denotes u"V v , and 

= V M it M is the volume expansion rate of the it M worldlines, 
A^ = Ufj, = Ai^s is their 4- acceleration, 

Cjui/ = Vuu v \ is their shear rate, 
oj^ = — icurlitju = cj/^n is their vorticity rate. 
On a Friedmann brane, 

^ = ^ = ^ = 0, 6 = 3if, (3.64) 

where = a/a is the Hubble rate. The covariant spatial curl is given by 

curl V M = e^p^V , curlW^ = e af n l $ a W fi v) , (3.65) 

where e^ a p is the projection orthogonal to of the 4D brane alternating tensor, and V M is the projected part of the 
brane covariant derivative, defined by 

V^-.-./j = (W^F a -... p ) ±u = h^h a 7 ■ ■ ■ hp 5 V v F^-... s . (3.66) 

In a local inertial frame at a point on the brane, with u M = S^o, we have: = Aq = ujq = cro^ = £oap = curlVo = 
curlWop and 

V^-F"-...^ = 6^5° j ■ ■■6 /3 k S/ i F j -... k (local inertial frame) , (3.67) 

where i, j, k = 1, 2, 3. 

The absence of bulk source terms in the conservation equations is a consequence of having A5 as the only 5D source 
in the bulk. For example, if there is a bulk scalar field, then there is energy-momentum exchange between the brane 
and bulk (in addition to the gravitational interaction) [29l |4C| . 

Equation l|3.30|l may be called the "nonlocal conservation equation" . Projecting along m m gives the nonlocal energy 
conservation equation, which is a propagation equation for pg. In the general, nonlinear case, this gives 

pe + \&Ps + V% + 2A*<% + a^n^ 
= [67^7^ + 6(p + pK%„ + 29 (2<z% + tt^tt^) + 2A"g%„ 

- 4<f V M p + ^V%„ + T^Vrfv - 2<j^ir afl Tr„ a - 2a^q^] . (3.68) 
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Projecting into the comoving rest space gives the nonlocal momentum conservation equation, which is a propagation 
equation for q^: 

= [-^(p + p)W^p + 6(p + P )W u TT fw + q u TT { ^ ) +7r M "V„(2p + 5p) 

2 / \ 28 

- ^k 010 (V^ap + 3V a ir Pll ) - 3 7 r AlQ V (3 7r Q ' 3 + —q u V^q v 

+ ipA"^ - Sn^A^ + ! Vr^TTa/s - ir^o^qp 

+ <y^ aP qp + Tr^e^Uaqp - e m pu) a 7r 0u q v + 4(p + p)<dq^ 

+ 6q ti A"q, + ^-A^q, + Aq^ a0 Tr a f 3 . (3.69) 

The 1+3-covariant decomposition shows two key features: 

• inhomogeneous and anisotropic effects from the 4D matter-radiation distribution on the brane are a source for 
the 5D Weyl tensor, which nonlocally "backreacts" on the brane via its projection 8^; 

• there are evolution equations for the dark radiative (nonlocal, Weyl) energy (ps) and momentum (q £ ) densities 



(carrying scalar and vector modes from bulk gravitons), but there is no evolution equation for the dark radiative 
anisotropic stress (7rf„) (carrying tensor, as well as scalar and vector, modes), which arises in both evolution 
equations. 

In particular cases, the Weyl anisotropic stress 7T^ V may drop out of the nonlocal conservation equations, i.e., when 

we can neglect a^n^, and A^it^. This is the case when we consider linearized perturbations about an 

FRW background (which remove the first and last of these terms) and further when we can neglect gradient terms 
on large scales (which removes the second term). This case is discussed in Sec. VI. But in general, and especially in 
astrophysical contexts, the tt^ u terms cannot be neglected. Even when we can neglect these terms, tt^ arises in the 
field equations on the brane. 

All of the matter source terms on the right of these two equations, except for the first term on the right of Eq. I|3.(j9|l . 
are imperfect fluid terms, and most of these terms are quadratic in the imperfect quantities q^ and 7r M!/ . For a single 
perfect fluid or scalar field, only the V M /o term on the right of Eq. i|3.69|) survives, but in realistic cosmological and 
astrophysical models, further terms will survive. For example, terms linear in tt^u will carry the photon quadrupole 
in cosmology or the shear viscous stress in stellar models. If there are two fluids (even if both fluids are perfect), 
then there will be a relative velocity generating a momentum density q^ = pv^, which will serve to source nonlocal 
effects. 

In general, the 4 independent equations in Eqs. I|3.68|l and (|3.69() constrain 4 of the 9 independent components of 
£ M „ on the brane. What is missing, is an evolution equation for tt^ v , which has up to 5 independent components. 
These 5 degrees of freedom correspond to the 5 polarizations of the 5D graviton. Thus in general, the projection of 
the 5-dimensional field equations onto the brane does not lead to a closed system, as expected, since there are bulk 
degrees of freedom whose impact on the brane cannot be predicted by brane observers. The KK anisotropic stress 
7r^„ encodes the nonlocality. 

In special cases the missing equation does not matter. For example, if irz v — by symmetry, as in the case of 
an FRW brane, then the evolution of is determined by Eqs. I|3.(j8|l and I|3.(j9|l . If the brane is stationary (with 
Killing vector parallel to u M ), then evolution equations are not needed for £ M „, although in general will still be 
undetermined. However, small perturbations of these special cases will immediately restore the problem of missing 
information. 

If the matter on the brane has a perfect-fluid or scalar-field energy- momentum tensor, the local conservation 
equations (|3.62l) and l|3.63|l reduce to 

p + e(p + p)=0, (3.70) 
V tlP +{p + p)A ll =0, (3.71) 



16 



while the nonlocal conservation equations l|3.68|l and l|3.69[) reduce to 



PS + \&Pe + + 2A»qZ + a^n^ = , (3.72) 

QU + \®4 + ~V M p £ + \pzK + V u rf v + A^% + o^ql - e/ a io v q £ a 
(P + P) 



A 

Equation (|3~73|) shows that [13 



-V^p. (3.73) 



• if B^v = and the brane energy-momentum tensor has perfect fluid form, then the density p must be homoge- 
neous, V M/ o = 0; 

• the converse does not hold, i.e., homogeneous density does not in general imply vanishing E^. 

A simple example of the latter point is the FRW case: Eq. (|3.73|) is trivially satisfied, while Eq. I|3.72|l becomes 

P£ + AH P£ = . (3.74) 

This equation has the dark radiation solution 

P£= P£ o(^y ■ (3.75) 

If = 0, then the field equations on the brane form a closed system. Thus for perfect fluid branes with 
homogeneous density and — 0, the brane field equations form a consistent closed system. However, this is 
unstable to perturbations, and there is also no guarantee that the resulting brane metric can be embedded in a 
regular bulk. 

It also follows as a corollary that inhomogcncous density requires nonzero S^: 

V„p^0 => £„ v ^0. (3.76) 

For example, stellar solutions on the brane necessarily have £ M „ ^ in the stellar interior if it is non-uniform. 
Perturbed FRW models on the brane also must have £ M „ ^ 0. Thus a nonzero £^v, and in particular a nonzero 'K piin 
is inevitable in realistic astrophysical and cosmological models. 



E. Propagation and constraint equations on the brane 



The propagation equations for the local and nonlocal energy density and momentum density are supplemented by 
further 1+3-covariant propagation and constraint equations for the kinematic quantities <d, A^, lu^, u^, and for the 
free gravitational field on the brane. The kinematic quantities govern the relative motion of neighbouring fundamental 
world-lines. The free gravitational field on the brane is given by the brane Weyl tensor G^ap. This splits into the 
gravito-electric and gravito-magnetic fields on the brane: 

E^u — C^avfivP'vP — , Hfj, v = -e ^ a pC a P v^u 1 — , (3.77) 

where E^ u is not to be confused with £^ v . The Ricci identity for 

VfrVv]U a = -Rc a ,tfU f> , (3-78) 

and the Bianchi identities 

VC^ - V [(U (-R v]a + \Rg v] ^j , (3.79) 

produce the fundamental evolution and constraint equations governing the above covariant quantities. The field equa- 
tions are incorporated via the algebraic replacement of the Ricci tensor R^ v by the effective total energy-momentum 
tensor, according to Eq. (|3.2U|) . The brane equations are derived directly from the standard general relativity versions 



17 



by simply replacing the energy- momentum tensor terms p, . . . by ptot, ■ ■ • ■ For a general fluid source, the equations 
are given m |3g. In the case of a single perfect fluid or minimally-coupled scalar field, the equations reduce to the 
following nonlinear equations. 

Generalized Raychaudhuri equation (expansion propagation): 

1 



Vorticity propagation: 



6 + -e 2 + - 2lu^ - VM M + A^A" + — (p + 3p) - A 



— (2p + 3p)^- K 2 P£ . 



2 1 



Shear propagation: 



TV- 



Gravito-electric propagation (Maxwell- Weyl E-dot equation): 



(3.80) 



(3.81) 



(3.82) 



+ ©-E/w ™ curl il^ + — (p + p)o>„ 
— 2A a e a p^H u - } 13 — 3a a ^E^ a + uj a e a p^E v ) 



+ 6A (A1 ^ + 3a Q <Al ^ )Q + 3w Q £ Q/3(Al7 r^ . 
Gravito-magnetic propagation (Maxwell- Weyl H-dot equation): 

H^ u) + QH^ + curl£ MW - 3cr a</1 Jf I , ) Q + w a e a/ j( M iT„) /5 + 2A a e aP{ll E„/ 

- 2 r i 
curl tt^ - Suj^q^ + a a (p,£v) a0 q$ ■ 



2 L 



Vorticity constraint: 



Shear constraint: 



Gravito-magnetic constraint: 



Vw„ - = . 



curl^ - - V^e + 2e uua uj v A a 



2 £ 



curlo^ + V( M Wj,) - i!^ + 2^4^^) = . 
Gravito-electric divergence (Maxwell- Weyl div-E equation) : 

2 2 

= y ^V^p + y (2V mP£ - 29^ - 3V^7T^ + 3a/ - 9e/ a w^ 
Gravito-magnetic divergence (Maxwell- Weyl div-H equation): 



izuj^ - 3curlg£ - 3e/ a a v l3 'K^ - 3?r^w 



(3.83) 



(3.84) 



(3.85) 



(3.86) 



(3.87) 



(3- 



(3.89) 



Gauss-Codazzi equations on the brane (with lu^ = 0): 



2 o 2 

R + o 02 - - 2k2 P - 2A = k 2 ^- + 2k 2 P£ , 

3 A 

where Rt v is the Ricci tensor for 3-surfaces orthogonal to u M on the brane and R = WR^. 
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FIG. 4: The evolution of the dimensionless shear parameter S7 s hcar = a 2 /6H 2 on a Bianchi I brane, for a V 
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model. 

and late-time expansion of the universe is isotropic, but the shear dominates during an intermediate anisotropic stage 



±m 2 < 



The standard 4D general relativity results are regained when A -1 — > and £ M „ = 0, which sets all right hand 
sides to zero in Eqs. (|3.80l) ~ (|3.91|) . Together with Eqs. H3.70fl - H3.73l) . these equations govern the dynamics of the 
matter and gravitational fields on the brane, incorporating both the local, high-energy (quadratic energy-momentum) 
and nonlocal, KK (projected 5D Weyl) effects from the bulk. High-energy terms are proportional to p/A, and are 
significant only when p > A. The KK terms contain p£, and 7rf„, with the latter two quantities introducing 
imperfect fluid effects, even when the matter has perfect fluid form. 

Bulk effects give rise to important new driving and source terms in the propagation and constraint equations. The 
vorticity propagation and constraint, and the gravito-magnetic constraint have no direct bulk effects, but all other 
equations do. High-energy and KK energy density terms are driving terms in the propagation of the expansion 0. 
The spatial gradients of these terms provide sources for the gravito-electric field E^ v . The KK anisotropic stress 
is a driving term in the propagation of shear a^ v and the gravito-electric/ -magnetic fields, E/H^, and the KK 
momentum density is a source for shear and the gravito-magnetic field. The 4D Maxwell- Weyl equations show in 
detail the contribution to the 4D gravito-electromagnetic field on the brane, i.e., (E^, H^), from the 5D Weyl field 
in the bulk. 

An interesting example of how high-energy effects can modify general relativistic dynamics arises in the analysis of 
isotropization of Bianchi spacetimes. For a Binachi type I brane, Eq. (|3.91|) becomes [3^| 



*-T'( 1 + s) + F- 



(3.92) 



if we neglect the dark radiation, where a and H are the average scale factor and expansion rate, and S is the shear 
constant. In general relativity, the shear term dominates as a — > 0, but in the brane- world, the high-energy p 2 term 
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will dominate if w > 0, so that the matter-dominated early universe is isotropic |38l l4ll |42| . This is illustrated in 
Fig. 4. 

Note that this conclusion is sensitive to the assumption that p£ ss 0, which, by Eq. H3.72[l implies the restriction 

<7""7r£„«0. (3.93) 

Relaxing this assumption can lead to non-isotropizing solutions |4flj . 

The system of propagation and constraint equations, i.e. Eqs. <|3.7L)I) - (|3.73|) and (|3.8U|) (|3.91l) . is exact and nonlinear, 
applicable to both cosmological and astrophysical modelling, including strong-gravity effects. In general the system 
of equations is not closed: there is no evolution equation for the KK anisotropic stress tt^. 



IV. GRAVITATIONAL COLLAPSE AND BLACK HOLES ON THE BRANE 



The physics of brane-world compact objects and gravitational collapse is complicated by a number of factors, 
especially the confinement of matter to the brane, while the gravitational field can access the extra dimension, and 
the nonlocal (from the brane viewpoint) gravitational interaction between the brane and the bulk. Extra-dimensional 
effects mean that the 4D matching conditions on the brane, i.e., continuity of the induced metric and extrinsic 
curvature across the 2-surface boundary, are much more complicated to implement |44l EEj . High-energy corrections 
increase the effective density and pressure of stellar and collapsing matter. In particular this means that the effective 
pressure does not in general vanish at the boundary 2-surface, changing the nature of the 4D matching conditions 
on the brane. The nonlocal KK effects further complicate the matching problem on the brane, since they in general 
contribute to the effective radial pressure at the boundary 2-surface. Gravitational collapse inevitably produces 
energies high enough, i.e., p^S> A, to make these corrections significant. 

We expect that extra-dimensional effects will be negligible outside the high-energy, short-range re gim e. The cor- 
rections to the weak-field potential, Eq. I|2.23|l . are at the second post-Newtonian (2PN) level [47|. However, 
modifications to Hawking radiation may bring significant corrections even for solar-sized black holes, as discussed 
below. 

A vacuum on the brane, outside a star or black hole, satisfies the brane field equations 

R la> = -£^, R\ = = £\, V^ = 0. (4.1) 

The Weyl term E^ v will carry an imprint of high-energy effects that source KK modes (as discussed above). This 
means that high-energy stars and the process of gravitational collapse will in general lead to deviations from the 4D 
general relativity problem. The weak-field limit for a static spherical source, Eq. I|2.23|l . shows that £ M „ must be 
nonzero, since this is the term responsible for the corrections to the Newtonian potential. 

A. The black string 

The projected Weyl term vanishes in the simplest candidate for a black hole solution. This is obtaine d by assuming 
the exact Schwarzschild form for the induced brane metric and "stacking" it into the extra dimension |48j 

= e-W/'fadafdar + dy 2 , (4.2) 

~9,u = e 2|y|/ V = -(1 - 2GM/r)dt 2 + 1 _ 2 + rW . (4.3) 

(Note that Eq. (|4.2|) is in fact a solution of the 5D field equations (|2.4|l if is any 4D Einstein vacuum solution, 
i.e., if R^v = 0, and this can be generalized to the case = —Ag^ [43,lfiCl|.) 

Each {y = const } surface is a 4D Schwarzschild spacetime, and there is a line singularity along r = for all y. 
This solution is known as the Schwarzschild black string, which is clearly not localized on the brane y — 0. Although 
^Cabcd 7^ 0, the projection of the bulk Weyl tensor along the brane is zero, since there is no correction to the 4D 
gravitational potential: 

V(r) = — => 8^ = 0. (4.4) 
r 

The violation of the perturbative corrections to the potential signals some kind of non-AdS5 pathology in the bulk. 
Indeed, the 5D curvature is unbounded at the Cauchy horizon, as y — > oo |48j : 

^Rabcd ^R ABCD = f + . (4.5) 
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Furthermore, the black string is unstable to large-scale perturbations [5l|. 

Thus the "obvious" approach to finding a brane black hole fails. An alternative approach is to seek solutions of 



the brane field equations with nonzero £ 



Brane solutions of static black hole exteriors with 5D corrections 



to the Schwarzschild metric have been found [U El E3- but the bulk metric for these solutions has not been 
found. Numerical integration into the bulk, starting from static black hole solutions on the brane, is plagued with 
difficulties IE3. 



B. Taylor expansion into the bulk 



One can use a Taylor expansion, as in Eq. (|3.39|) , in order to probe properties of a static black hole on the brane \54 
(An alternative expansion scheme is discussed in |55j.) For a vacuum brane metric, 



9tiv{x,y) = g^(x,0) - £ t _ tu (x,0+)y 2 - j£^{x,0+)\y\ 3 
1 r 32 

~t~Y^ ^—^£^iu ~^2^'t J:l/ ~t~ Q£[t £qlv 



y=o+ 



y 



(4.6) 



This shows in particular that the propagating effect of 5D gravity arises only at the fourth order of the expansion. 
For a static spherical metric on the brane, 



dr 2 



g^dx^dx" = -F(r)dt 2 + + r 2 dQ 2 



H(r) 



the projected Weyl term on the brane is given by 

F r , 
£00 = — H — 



1 -H 



f — 

rtl 



F' 1 - H 
~¥ r 



r /F' H' 
See - -^H + -H(- + - 



(4.7) 

(4.8) 
(4.9) 
(4.10) 



These components allow one to evaluate the metric coefficients in Eq. I|4.6() . For example, the area of the 5D horizon 
is determined by gee] defining ip(r) as the deviation from a Schwarzschild form for H, i.e., 



/ n 2m , . . 

H(r) = 1 + V(r) , 

r 



where m is constant, we find 



!l,,n(r. >i> =r 2 [l + j\y\\y 2 + 



(4.11) 



(4.12) 



This shows how -0 and its r-derivatives determine the change in area of the horizon along the extra dimension. For 
the black string, ip = and we have ggg(r,y) — r 2 . For a large black hole, with horizon scale ^> I, we have from 
Eq. lt2~23|l that 



Ami 2 
' 3r 3 



(4.13) 



This implies that gee is decreasing as we move off the brane, consistent with a pancake-like shape of the horizon. 
However, note that the horizon shape is tubular in Gaussian normal coordinates j56j. 

C. The "tidal charge" black hole 

The equations 14.1(1 form a system of constraints on the brane in the stationary case, including the static spherical 
case, for which 



6 = = U)^ = , p E = = q £ = tt 



(4.14) 



21 



The nonlocal conservation equations VE^v = reduce to 

1^ 4 



3 p " 3 

where, by symmetry, 



V M p £ + -P£^ + V v k% + A v -k% = , (4.15) 



= n £ ( i v - r ^ ) . ( 4 - 16 ) 



for some n,c-(r), with r M the unit radial vector. The solution of the brane field equations requires the input of E^ v 

from the 5D solution. In the absence of a 5D solution, one can make an assumption about £ ^ or to close the 4D 
equations . 

If we assume a metric on the brane of Schwarzschild-like form, i.e., H = F in Eq. I|4.7J) . then the general solution 
of the brane field equations is 

F = 1 _2^M + 2^Q ) (41?) 
2G£Q 

£^u = ^j— [u^u v - Ir^rv + h^v] , (4-18) 

where Q is a constant. It follows that the KK energy density and anisotropic stress scalar (defined via Eq. (|4.16|l ^ are 
given by 

P e = ^- i =\n £ . (4.19) 
47r r 4 2 

The solution Eq. (|4.17(l has the form of the general relativity Reissner-Nordstrom solution, but there is no electric 
field on the brane. Instead, the nonlocal Coulomb effects imprinted by the bulk Weyl tensor have induced a "tidal" 
charge parameter Q, where Q — Q(M), since M is the source of the bulk Weyl field. We can think of the gravitational 
field of M being "reflected back" on the brane by the negative bulk cosmological constant [57j • If we impose the 
small-scale perturbative limit (r <C £) in Eq. (|2.22|) . we find that 

Q = -2M. (4.20) 

Negative Q is in accord with the intuitive idea that the tidal charge strengthens the gravitational field, since it arises 
from the source mass M on the brane. By contrast, in the Reissner-Nordstrom solution of general relativity, Q ot +q 2 , 
where q is the electric charge, and this weakens the gravitational field. Negative tidal charge also preserves the spacelike 
nature of the singularity, and it means that there is only one horizon on the brane, outside the Schwarzschild horizon: 



r h = GM 



1 - 



21Q 
GM 2 



= GM 



41 
GM 



(4.21) 



The tidal-charge black hole metric does not satisfy the far-field r~ 3 correction to the gravitational potential, as 
in Eq. I|2.23|l . and therefore cannot describe the end-state of collapse. However, Eq. (|4.17|l shows the correct 5D 
behaviour of the potential (oc r~ 2 ) at short distances, so that the tidal-charge metric could be a good approximation 
in the strong-field regime for small black holes. 



D. Realistic black holes 



Thus a simple brane-based approach, while giving useful insights, does not lead to a realistic black hole solution. 
There is no known solution representing a realistic black hole localized on the brane, which is stable and without 
naked singularity. This remains a key open question of nonlinear brane-world gravity. (Note that an exact solution is 
known for a black hole on a 1+2-brane in a 4D bulk but this is a very special case.) Given the nonlocal nature 
of £ M „, it is possible that the process of gravitational collapse itself leaves a signature in the black hole end-state, in 
contrast with general relativity and its no- hair theorems. There are contradictory indications about the nature of the 
realistic black hole solution on the brane: 
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• numerical simulations of highly relativistic static stars on the brane [59] indicate that general relativity remains 
a good approximation; 

• exact analysis of Oppenheimer-Snyder collapse on the brane shows that the exterior is non-static |45|. and this 
is extended to general collapse by arguments based on a generalized AdS/CFT correspondence |6(ll6l|. 

The first result suggests that static black holes could exist as limits of increasingly compact static stars, but the 
second result and conjecture suggest otherwise. This remains an open question. More recent numerical evidence is 
also not conclusive, and it introduces further possible subtleties to do with the size of the black hole 

On very small scales relative to the AdSs curvature scale, r I, the gravitational potential becomes 5D, as shown 
in Eq. 

VW-S£-3£. (4-22) 

In this regime, the black hole is so small that it does not "see" the brane, so that it is approximately a 5D Schwarzschild 
(static) solution. However, this is always an approximation because of the self-gravity of the brane (the situation is 
different in ADD-type brane- worlds where there is no brane tension). As the black hole size increases, the approxima- 
tion breaks down. Nevertheless, one might expect that static solutions exist on sufficiently small scales. Numerical 
investigations appear to confirm this |62|: static metrics satisfying the asymptotic AdSs boundary conditions are 
found if the horizon is small compared to £, but no numerical convergence can be achieved close to I. The numerical 
instability that sets in may mask the fact that even the very small black holes are not strictly static. Or it may be 
that there is a transition from static to non-static behaviour. Or it may be that static black holes do exist on all 
scales. 

The 4D Schwarzschild metric cannot describe the final state of collapse, since it cannot incorporate the 5D behaviour 
of the gravitational potential in the strong-field regime (the metric is incompatible with massive KK modes). A non- 
perturbative exterior solution should have nonzero E^u in order to be compatible with massive KK modes in the 
strong-field regime. In the end-state of collapse, we expect a E^ v which goes to zero at large distances, recovering the 
Schwarzschild weak-field limit, but which grows at short range. Furthermore, may carry a Weyl "fossil record" 
of the collapse process. 

E. Oppenheimer-Snyder collapse gives a non-static black hole 

The simplest scenario in which to analyze gravitational collapse is the Oppenheimer-Snyder model, i.e., collapsing 
homogeneous and isotropic dust |45j . The collapse region on the brane has FRW metric, while the exterior vacuum 
has an unknown metric. In 4D general relativity, the exterior is a Schwarzschild spacetime; the dynamics of collapse 
leaves no imprint on the exterior. 



The collapse region has metric 



, a( T ) 2 \dr 2 +r 2 ain 2 } , 
ds 2 = -dr 2 + K -L± - g— -„ J , 4.23 
(1 + kr 2 /4) 2 y ' 



where the scale factor satisfies the modified Friedmann equation (see below) 

a 2 8irG f p ps 



P 1 + + — • ( 4 - 24 ) 



a 2 3 r \~ ' 2A ' p 
The dust matter and the dark radiation evolve as 

P = Po(j) 3 . Ps = Pso{^)\ (4.25) 

where ao is the epoch when the cloud started to collapse. The proper radius from the centre of the cloud is R(t) = 
ra(r)/(l + jkr 2 ). The collapsing boundary surface £ is given in the interior comoving coordinates as a free-fall 
surface, i.e. r = r = const, so that i?s(r) = roa(r)/(l + \kr^). 

We can rewrite the modified Friedmann equation on the interior side of £ as 

• a _ 2GM 3GM 2 Q 

R - — + 4^ + R 2 - +E > (4 - 26) 
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where the "physical mass" M (total energy per proper star volume), the total "tidal charge" Q and the "energy" per 
unit mass are given by 

- 3(1 + Ifcr-S) 3 ' 1 ' 

Q ~ oW ,4 ' 28) 

Now we assume that the exterior is static, and satisfies the standard 4D junction conditions. Then we check whether 
this exterior is physical by imposing the modified Einstein equations (14.11) . We will find a contradiction. 

The standard 4D Darmois-Israel matching conditions, which we assume hold on the brane, require that the metric 
and the extrinsic curvature of £ be continuous (there are no intrinsic stresses on £). The extrinsic curvature is 
continuous if the metric is continuous and if R is continuous. We therefore need to match the metrics and R across 
S. 

The most general static spherical metric that could match the interior metric on £ is 



ds 2 = -F{R) 2 



2Gm{R) 
1 R 



dt2+ l-2GMR)/R +R2dn2 - (4 ' 30) 



We need two conditions to determine the functions F(R) and m{R). Now £ is a freely falling surface in both 
metrics, and the radial geodesic equation for the exterior metric gives R 2 = — 1 + 2Gm(R) / R + E/F(R) 2 , where 
E is a constant and the dot denotes a proper time derivative, as above. Comparing this with Eq. I|4.26|l gives one 
condition. The second condition is easier to derive if we change to null coordinates. The exterior static metric, with 
dv = dt + dR/[(l - 2Gm/R)F], becomes 

ds 2 = -F 2 (l-2Gm/R)dv 2 + 2FdvdR + R 2 dtt 2 . (4.31) 

The interior Robertson- Walker metric takes the form 

„ t\ \\ - (k + a 2 )R 2 /a 2 } dv 2 2t v dvdR , , , 

ds 2 = — ^ ji , ' . ' J + - '= =+R 2 dn 2 , 4.32 

(l~kR 2 /a 2 ) ^l-kR 2 /a 2 V ' 

where dr — T tV dv + (1 + \kr 2 )dR/[ra — 1 + jkr 2 ]. Comparing Eqs. (|4.31|) and 14.32fl on E gives the second condition. 
The two conditions together imply that F is a constant, which we can take as F(R) = 1 without loss of generality 
(choosing E = E + 1), and that 

3/W 2 (1 

m ^- M+ S^ + 2^R- ^ 

In the limit A -1 = = Q, we recover the 4D Schwarzschild solution. In the general brane- world case, Eqs. (|4.30(l and 
(|4.33|) imply that the brane Ricci scalar is 

R> (4 34) 

However, this contradicts the field equations (|4.1|l . which require 

R% = 0. (4.35) 

It follows that a static exterior is only possible if Mj A = , which is the 4D general relativity limit. In the brane- world, 
collapsing homogeneous and isotropic dust leads to a non-static exterior. Note that this no-go result does not require 
any assumptions on the nature of the bulk spacetime, which remains to be determined. 

Although the exterior metric is not determined (see [63j for a toy model), we know that its non-static nature arises 
from 



• 5D bulk graviton stresses, which transmit effects nonlocally from the interior to the exterior, and 
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• the non-vanishing of the effective pressure at the boundary, which means that dynamical information from the 
interior can be conveyed outside via the 4D matching conditions. 

The result suggests that gravitational collapse on the brane may leave a signature in the exterior, dependent upon 
the dynamics of collapse, so that astrophysical black holes on the brane may in principle have KK "hair" . It is possible 
that the non-static exterior will be transient, and will tend to a static geometry at late times, close to Schwarzschild 
at large distances. 



F. AdS/CFT and black holes on 1-brane RS-type models 

Oppenheimer-Snyder collapse is very special; in particular, it is homogeneous. One could argue that the non-static 
exterior arises because of the special nature of this model. However, the underlying reasons for non-static behaviour 
are not special to this model; on the contrary, the role of high-energy corrections and KK stresses will if anything 
be enhanced in a general, inhomogeneous collapse. There is in fact independent heuristic support for this possibility, 
arising from the AdS/CFT correspondence. 

The basic idea of the correspondence is that the classical dynamics of the AdSs gravitational field correspond to 
the quantum dynamics of a 4D conformal field theory on the brane. This correspondence holds at linear perturbative 
order so that the RS 1-brane infinite AdS$ brane- world (without matter fields on the brane) is equivalently 
described by 4D general relativity coupled to conformal fields, 

<V = SnGT^ . (4.36) 

According to a conjecture the correspondence holds also in the case where there is strong gravity on the brane, 
so that the classical dynamics of the bulk gravitational field of the brane black hole are equivalent to the dynamics of 
a quantum-corrected 4D black hole (in the dual CFT-plus-gravity description) . In other words E3 j 



quantum backreaction due to Hawking radiation in the 4D picture is described as classical dynamics in the 5D 
picture; 



• the black hole evaporates as a classical process in the 5D picture, and there is thus no stationary black hole 
solution in RS 1-brane. 

A further remarkable consequence of this conjecture is that Hawking evaporation is dramatically enhanced, because 
of the very large number, of order (£/£ p ) 2 , of CFT modes. The energy loss rate due to evaporation is 

M ( 1 \ , 

where N is the number of light degrees of freedom. Using N ~ £ 2 /G, this gives an evaporation timescale [6(j 

M \ 3 / 1 mm \ 2 



[—) (4 - 38) 

A more detailed analysis [(54| shows that this expression should be multiplied by a factor ss 100. Then the existence 
of stellar-mass black holes on long time scales places limits on the AdSs curvature scale that are more stringent than 
the table-top limit, Eq. I|1.6fl . The existence of black hole X-ray binaries implies 

f<10~ 2 mm, (4.39) 

already an order of magnitude improvement on the table-top limit. 

One can also relate the Oppenheimer-Snyder result to these considerations. In the AdS/CFT picture, the non- 
vanishing of the Ricci scalar, Eq. I|4.34|) arises from the trace of the Hawking CFT energy-momentum tensor, as in 
Eq. 1)4. 36[) . If we evaluate the Ricci scalar at the black hole horizon, R ~ 2GM, using A = 6M®/M 2 , we find 



mi 2 e 

M 4 

The CFT trace on the other hand is given by ~ NT^/M 2 so that 



(4.41) 
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Thus the Oppenheimer-Snyder result is qualitatively consistent with the AdS/CFT picture. 

Clearly the black hole solution, and the collapse process that leads to it, have a far richer structure in the brane- world 
than in general relativity, and deserve further attention. In particular, two further topics are of interest: 

• Primordial black holes in 1-brane RS-type cosmology have been investigated in 0, Ell ■ High-energy effects 
in the early universe (see the next section) can significantly modify the evaporation and accretion processes, 
leading to a prolonged survival of these black holes. Such black holes evade the enhanced Hawking evaporation 
described above when they are formed, because they are much smaller than t. 

• Black holes will also be produced in particle collisions at energies > M5, possibly well below the Planck scale. 
In ADD brane- worlds, where M^+d — O(TeV) is not ruled out by current observations if d > 1, this raises the 
exciting prospect of observing black hole production signatures in the next-generation colliders and cosmic ray 
detectors (see 0,|6^). 

V. BRANE- WORLD COSMOLOGY: DYNAMICS 

A 1+4-dimensional spacetime with spatial 4-isotropy (4D spherical/ plane/ hyperbolic symmetry) has a natural 
foliation into the symmetry group orbits, which are 1+3-dimensional surfaces with 3-isotropy and 3-homogeneity, i.e. 
FRW surfaces. In particular, the AdSs bulk of the RS brane-world, which admits a foliation into Minkowski surfaces, 
also admits an FRW foliation since it is 4-isotropic. Indeed this feature of 1-brane RS-type cosmological brane-worlds 
underlies the importance of the AdS/CFT correspondence in brane-world cosmology |l3j . 

The generalization of AdSs that preserves 4-isotropy and solves the vacuum 5D Einstein equation 1)2. 4J) . is 
Schwarzschild-AdSs, and this bulk therefore admits an FRW foliation. It follows that an FRW brane-world, the 
cosmological generalization of the RS brane-world, is a part of Schwarzschild-AdSs, with the ^-symmetric FRW 
brane at the boundary. (Note that FRW branes can also be embedded in non-vacuum generalizations, e.g. in 
Reissner-Nordstrom-AdSs and Vaidya-AdSs.) 

In natural static coordinates, the bulk metric is 



^ds 2 



-F(R)dT 



F(R) = K 



e 2 



dR 2 
F(R) 



R 



dr 2 



1 - Kr 2 



l dVL' 



m 
R? 



(5.1) 
(5.2) 



where K = 0, ±1 is the FRW curvature index and m is the mass parameter of the black hole at R = (recall that 
the 5D gravitational potential has R~ 2 behaviour). The bulk black hole gives rise to dark radiation on the brane via 
its Coulomb effect. The FRW brane moves radially along the 5th dimension, with R = a(T), where a is the FRW 
scale factor, and the junction conditions determine the velocity via the Friedmann equation for a 36]. Thus one can 
interpret the expansion of the universe as motion of the brane through the static bulk. In the special case m = and 
da/dT — 0, the brane is fixed and has Minkowski geometry, i.e., the original RS 1-brane brane-world is recovered, in 
different coordinates. 

The velocity of the brane is coordinate-dependent, and can be set to zero. We can use Gaussian normal coordinates, 
in which the brane is fixed but the bulk metric is not manifestly static [24|: 



(5) ds 2 = -N 2 {t,y)dt 2 + A 2 (t,y) 



dr 2 



1 - Kr 2 



dy 2 



(5.3) 



Here a(t) = A(t, 0) is the scale factor on the FRW brane at y = 0, and t may be chosen as proper time on the brane, 
so that N(t, 0) = 1. In the case where there is no bulk black hole (m = 0), the metric functions are 



N 
A 



M*,y) 

a(t) 



o(t) 



cosh 



(?H^Mt) 



Again, the junction conditions determine the Friedmann equation. The extrinsic curvature at the brane is 

„ M ,. ( N' A' A' A' 



(5.4) 
(5.5) 

(5.6) 
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Then, by Eq. (|3~T§|) . 



T 

The field equations yield the first integral [24| 



(2p + 3p-A), (5.7) 

brane 

= -f(P + A). (5.8) 

branc D 



(AA') 2 -| I i 2 + ^ 4 + m = 0, (5.9) 

where m is constant. Evaluating this at the brane, using Eq. (|5.8|l . gives the modified Friedmann equation l|5.11[) . 

The dark radiation carries the imprint on the brane of the bulk gravitational field. Thus we expect that £ M „ for the 
Friedmann brane contains bulk metric terms evaluated at the brane. In Gaussian normal coordinates (using the field 
equations to simplify), 



A" 

£0 °=4 



4p r, = -(ifV) (-in, 



Either form of the cosmological metric, Eq. (|5.1() or l|5.3(l . may be used to show that 5D gravitational wave signals 
can take "short-cuts" through the bulk in travelling between points A and B on the brane [g^. The travel time for 
such a graviton signal is less than the time taken for a photon signal (which is stuck to the brane) from A to B. 

Instead of using the junction conditions, we can use the covariant 3D Gauss-Codazzi equation l|3.91|l to find the 
modified Friedmann equation: 



H =TK 1+ 2T) + ^ + 3 A -^ (5 ' n) 



on using Eq. (|3.75|l . where 



m=^—p £0 al. (5.12) 



The covariant Raychauhuri equation l|3.80|l yields 



2 yr 1 ' \ A 



a 



which also follows from differentiating Eq. i|5.11[l and using the energy conservation equation. 

When the bulk black hole mass vanishes, the bulk geometry reduces to AdSs and p£ = 0. In order to avoid a naked 
singularity, we assume that the black hole mass is non-negative, so that pg > 0. [By Eq. 1)5. 2|) . it is possible to avoid 
a naked singularity with negative m when K = —1, provided \m\ < £ 2 /4.] This additional effective relativistic degree 
of freedom is constrained by nucleosynthesis and CMB observations to be no more than ~5% of the radiation energy 
density [HHHz^ 



Ps 

Prad 



<0.05 (5.14) 



The other modification to the Hubble rate is via the high-energy correction p/X. In order to recover the observational 
successes of general relativity, the high-energy regime where significant deviations occur must take place before 
nucleosynthesis, i.e., cosmological observations impose the lower limit 

A > (1 MeV) 4 => M 5 > 10 4 GeV. (5.15) 

This is much weaker than the limit imposed by table-top experiments, Eq. 12.24JI . Since p 2 /X decays as a -8 during 
the radiation era, it will rapidly become negligible after the end of the high-energy regime, p = A. 

If ps = and K = = A, then the exact solution of the Friedmann equations for w = p/ p = const is |24| 

a= const [t(i + t A )] 1 / 3(w+1 ), t\ — ¥ p < (l + io^lfj^sec, (5.16) 
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where w > —I. If ps (but K = = A), then the solution for the radiation era (w = |) is [6jJ 

a = const [t(t + tA)] 1 / 4 , t x = . (5.17) 

4V7tA(1 + p e /p) 

For t t\ we recover from Eqs. (|5.16(1 and l|5.17(l the standard behaviour, a cx t 2 / 3 ( w+1 > , whereas for i <C t\, we have 
the very different behaviour of the high-energy regime, 

When io = —1 we have p = po from the conservation equation. If K = = A, we recover the de Sitter solution for 
ps = and an asymptotically de Sitter solution for ps > 0: 



p s = 0: o = aoexp[£r (<-*o)], #o - «\/ ( 1 + £f ) - ( 5 - 19 ) 
/ Of>0: a 2 = y^sinh[2ff (*-*o)]- (5-20) 



A qualitative analysis of the Friedmann equations is given in [41 



A. Brane-world inflation 



In 1-brane RS-type brane-worlds, where the bulk has only a vacuum energy inflation on the brane must be driven 
by a 4D scalar field trapped on the brane. In more general brane-worlds, where the bulk contains a 5D scalar field, 
it is possible that the 5D field induces inflation on the brane via its effective projection |7_l|. 

More exotic possibilities arise from the interaction between two branes, including possible collision, which is mediated 
by a 5D scalar field and which can induce either inflation |72| or a hot big-bang radiation era, as in the "ekpyrotic" or 
cyclic scenario 73], or in colliding bubble scenarios 

[Z3- (See also HI for collidingb ranes in an M theory approach.) 
Here we discuss the simplest case of a 4D scalar field <f> with potential V(<f>) (see |7(| for a review). 

High-energy brane-world modifications to the dynamics of inflation on the brane have been investigated [38l FtI l78| . 
Essentially the high-energy corrections provide increased Hubble damping, since p> A implies H is larger for a given 
energy than in 4D general relativity. This makes slow-roll inflation possible even for potentials that would be too 
steep in standard cosmology [23, [23, • 

The field satisfies the Klein-Gordon equation 



4> + 3H<j> + V'{(f>) = 0. 



(5.21) 



In 4D general relativity, the condition for inflation, d > 0, is <j> 2 < V(tf>), i.e., p < — hp, where p = \<\> 2 + V and 
p = ^(f> 2 — V. The modified Friedmann equation leads to a stronger condition for inflation: using Eq. (|5.11() . with 
m = = A = K, and Eq. (|5.21|) . we find that 



a > 



w < 



1 + 2p/A 



1 + p/X 



(5.22) 



where the square brackets enclose the brane correction to the general relativity result. As pj\j — > 0, the 4D result 



w 



< — i is recovered, but for p > A, w must be more negative for inflation. In the very high-energy limit p/X 



oo, 



we have w < — |. When the only matter in the universe is a self- interacting scalar field, the condition for inflation 
becomes 



-V- 



V (h 



< 0, 



(5.23) 



which reduces to (j) 2 < V when p 
In the the slow-roll approximation, 



|0 2 + V < A. 



H 



K 2 

-XL 

3H 



(5.24) 
(5.25) 
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The brane-world correction term V/X in Eq. I|5.24() serves to enhance the Hubble rate for a given potential energy, 
relative to general relativity. Thus there is enhanced Hubble 'friction' in Eq. (|5.25|l . and brane-world effects will 
reinforce slow-roll at the same potential energy. We can see this by defining slow-roll parameters that reduce to the 
standard parameters in the low-energy limit: 
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(5.26) 
(5.27) 
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Self-consistency of the slow-roll approximation then requires e, \rj\ <C 1. At low energies, V <C A, the slow-roll 
parameters reduce to the standard form. However at high energies, V ^> A, the extra contribution to the Hubble 
expansion helps damp the rolling of the scalar field and the new factors in square brackets become ~ X/V: 

(5.28) 

where e gr , rj gv are the standard general relativity slow-roll parameters. In particular, this means that steep potentials 
which do not give inflation in general relativity, can inflate the brane-world at high energy and then naturally stop 
inflating when V drops below A. These models can be constrained because they typically end inflation in a kinetic- 
dominated regime and thus generate a blue spectrum of gravitational waves, which can disturb nucleosynthesis [791. 
They also allow for the novel possibility that the inflaton could act as dark matter or quintessence at low energies [73, 

m 

The number of e-folds during inflation, N = J Hdt, is, in the slow-roll approximation, 
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(5.29) 



Brane-world effects at high energies increase the Hubble rate by a factor V/2X, yielding more inflation between any 
two values of <f> for a given potential. Thus we can obtain a given number of e-folds for a smaller initial inflaton value 
4>i. For V > A, Eq. (|5~2^|l becomes 
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(5.30) 



The key test of any modified gravity theory during inflation, will be the spectrum of perturbations produced 
due to quantum fluctuations of the fields about their homogeneous background values. We will discuss brane-world 
cosmological perturbations in the next section. In general, perturbations on the brane are coupled to bulk metric 
perturbations, and the problem is very complicated. However on large scales on the brane, the density perturbations 
decouple from the bulk metric perturbations [35L l68l 18^ . For 1-brane RS-type models, there is no scalar zero- mode of 
the bulk graviton, and in the extreme slow-roll (de Sitter) limit, the massive scalar modes are heavy and stay in their 
vacuum state during inflation [8i| ■ Thus it seems a reasonable approximation in slow-roll to neglect the KK effects 
carried by when computing the density perturbations. 

To quantify the amplitude of scalar (density) perturbations we evaluate the usual gauge-invariant quantity 



P 



(5.31) 



which reduces to the curvature perturbation, 1Z, on uniform density hypersurfaces (5p = 0). This is conserved on 
large scales for purely adiabatic perturbations, as a consequence of energy conservation (independently of the field 
equations) |8^. The curvature perturbation on uniform density hypersurfaces is given in terms of the scalar field 
fluctuations on spatially flat hypersurfaces, Scj), by 



(5.32) 



The field fluctuations at Hubble crossing (k = aH) in the slow-roll limit are given by (Scj) 2 ) w (H/2n) , a result 
for a massless field in de Sitter space that is also independent of the gravity theory J84J. For a single scalar field 
the perturbations are adiabatic and hence the curvature perturbation £ can be related to the density perturbations 
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FIG. 5: The relation between the inflaton mass iti/Ma {Ma = M p ) and the brane tension (A/M4) 1 / 4 necessary to satisfy the 
COBE constraints. The straight line is the approximation used in Eq. 15.37L which at high energies is in excellent agreement 
with the exact solution, evaluated numerically in slow-roll. (From |77|.'l 



when modes re-enter the Hubble scale during the matter dominated era which is given by A 2 
slow- roll equations and Eq. H5.32[) . this gives 
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k=aH 



Thus the amplitude of scalar perturbations is increased relative to the standard result at a fixed value of <j> for a given 
potential. 

The scale-dependence of the perturbations is described by the spectral tilt 
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(5.34) 



where the slow-roll parameters are given in Eqs. H5.26[l and Ij5.27|l . Because these slow-roll parameters are both 
suppressed by an extra factor X/V at high energies, we see that the spectral index is driven towards the Harrison- 
Zel'dovich spectrum, n s — > 1, as V/\ — > 00; however, as explained below, this does not necessarily mean that the 
brane-world case is closer to scale-invariance than the general relativity case. 

As an example, consider the simplest chaotic inflation model V = ^m 2 (f) 2 . Equation l|5.29|l gives the integrated 
expansion from 0; to <p{ as 
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(5.35) 



The new high-energy term on the right leads to more inflation for a given initial inflaton value 4>i- 

The standard chaotic inflation scenario requires an inflaton mass m ~ 10 13 GeV to match the observed level of 
anisotropics in the cosmic microwave background (see below). This corresponds to an energy scale ~ 10 16 GeV when 
the relevant scales left the Hubble scale during inflation, and also to an inflaton field value of order 3M p . Chaotic 
inflation has been criticised for requiring super-Planckian field values, since these can lead to nonlinear quantum 
corrections in the potential. 

If the brane tension A is much below 10 16 GeV, corresponding to M5 < 10 17 GeV, then the terms quadratic in the 
energy density dominate the modified Friedmann equation. In particular the condition for the end of inflation given 
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in Eq. i|5.23[l becomes if 2 < |V. In the slow-roll approximation [using Eqs. I|5.24[l and l|5.25[) ] <j) ~ —M^/2n(j) and this 
yields 



In order to estimate the value of <$> when scales corresponding to large-angle anisotropics on the microwave background 
sky left the Hubble scale during inflation, we take N co \^ c ss 55 in Eq. I|5.35|l and <fi{ = <fi cn d- The second term on the 
right of Eq. (|5.35|) dominates, and we obtain 
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Imposing the COBE normalization on the curvature perturbations given by Eq. (|5.33(l requires 



in 2 \ m 4 ^ 5 



4.* ^ ^f^2xHT 5 . (5.38) 



45 / Mf 

Substituting in the value of </> co be given by Eq. i|5.37[l shows that in the limit of strong brane corrections, observations 
require 

to w 5 x 10~ 5 M 5 , cobc w 3 x 10 2 M 5 . (5.39) 

Thus for M5 < 10 17 GeV, chaotic inflation can occur for field values below the 4D Planck scale, 4> co bc < M p , although 
still above the 5D scale M 5 . The relation determined by COBE constraints for arbitrary brane tension is shown in 
Fig. 5, together with the high-energy approximation used above, which provides an excellent fit at low brane tension 
relative to M4. 

It must be emphasized that in comparing the high-energy brane-world case to the standard 4D case, we implicitly 
require the same potential energy. However, precisely because of the high-energy effects, large-scale perturbations 
will be generated at different values of V than in the standard case, specifically at lower values of V, closer to the 
reheating minimum. Thus there are two competing effects, and it turns out that the shape of the potential determines 
which is the dominant effect j8^- For the quadratic potential, the lower location on V dominates, and the spectral 
tilt is slightly further from scale invariance than in the standard case. The same holds for the quartic potential. Data 
from WMAP and 2dF can be used to constrain inflationary models via their deviation from scale invariance, and the 
high-energy brane-world versions of the quadratic and quartic potentials are thus under more pressure from data than 
their standard counterparts [83, as shown in Fig. 6. 

Other perturbation modes have also been investigated: 

• High-energy inflation on the brane also generates a zero-mode (4D graviton mode) of tensor perturbations, and 
stretches it to super-Hubble scales, as will be discussed below. This zero-mode has the same qualitative features 
as in general relativity, remaining frozen at constant amplitude while beyond the Hubble horizon. Its amplitude 
is enhanced at high energies, although the enhancement is much less than for scalar perturbations |8t| : 
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(5.41) 



Equation 15.41f) means that brane-world effects suppress the large-scale tensor contribution to CMB anisotropies. 
The tensor spectral index at high energy has a smaller magnitude than in general relativity, 

n t = -3e, (5.42) 

but remarkably the same consistency relation as in general relativity holds 80] : 

n t = -2^| . (5.43) 

This consistency relation persists when Z2 symmetry is dropped [s?! ] (and in a two-brane model with stabilized 
radion |s3)- ^ holds only to lowest order in slow-roll, as in general relativity, but the reason for this j8(| and 
the nature of the corrections 193 are not settled. 
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FIG. 6: Constraints from WMAP data on inflation models with quadratic and quartic potentials, where R is the ratio of tensor 
to scalar amplitudes and n is the scalar spectral index. The high energy (H.E.) and low energy (L.E.) limits are shown, with 
intermediate energies in between, and the 1-a and 2-a contours are also shown. (From [85f.1 



The massive KK modes of tensor perturbations remain in the vacuum state during slow-roll inflation mill. 
The evolution of the super-Hubble zero mode is the same as in general relativity, so that high-energy brane-world 
effects in the early universe serve only to rescale the amplitude. However, when the zero mode re-enters the 
Hubble horizon, massive KK modes can be excited. 

• Vector perturbations in the bulk metric can support vector metric perturbations on the brane, even in the 
absence of matter perturbations (see the next section). However, there is no normalizable zero mode, and the 
massive KK modes stay in the vacuum state during brane-world inflation [i^. Therefore, as in general relativity, 
we can neglect vector perturbations in inflationary cosmology. 

Brane-world effects on large-scale isocurvature perturbations in 2-field inflation have also been considered [93| . 
Brane-world (p)reheating after inflation is discussed in |94| . 



B. Brane-world instanton 



The creation of an inflating brane-world can be modelled as a de Sitter instanton in a way that closely follows the 
4D instanton, as shown in [26j. The instanton consists of two identical patches of AdSs joined together along a dS4 
brane with compact spatial sections. The instanton describes the "birth" of both the inflating brane and the bulk 
spacetime, which are together "created from nothing" , i.e. the point at the south pole of the de Sitter 4-sphere. The 
Euclidean AdS$ metric is 
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(5.44) 



where dfl 2 3 ^ is a 3-sphere, and r < rg. The Euclidean instanton interpolates between r = ("nothing") and r = tq 
(the created universe), which is a spherical brane of radius 



Hq 1 = £sinh(r /£) . 

After creation, the brane-world evolves according to the Lorentzian continuation, x ~ * iH^t + 7r/2, 



^ds 2 = dr 2 + (£H ) 2 smh 2 (r/£) 
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(5.45) 
(5.46) 
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FIG. 7: Brane-world instariton. (From |2£ 



Fig. 7.) 



C. Models with non-empty bulk 



le single-brane cosmological model can be generalized to include stresses other than A5 in the bulk. 

The simplest example arises from considering a charged bulk black hole, leading to the Reissner-Nordstrbm 
AdS 5 bulk metric [95|. This has the form of Eq. Q5.1JI . with 



(5.47) 



where q is the "electric" charge parameter of the bulk black hole. The metric is a solution of the 5D Einstein- 
Maxwell equations, so that ^Tab in Eq. 1)3.1)1 is the energy-momentum tensor of a radial static 5D "electric" 
field. In order for the field lines to terminate on the boundary brane, the brane should carry a charge —q. 
Since the RNAdSs metric is 4-isotropic, it is still possible to embed a FRW brane in it, which is moving in the 
coordinates of Eq. ()5.1)l . 

The effect of the black hole charge on the brane arises via the junction conditions and leads to the modified 
Friedmann equation |9fij . 
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(5.48) 



The field lines that terminate on the brane imprint on the brane an effective negative energy density —3q 2 / (k 2 o 6 ), 
which redshifts like stiff matter (w = 1). The negativity of this term introduces the possibility that at high 
energies it can bring the expansion rate to zero and cause a turn-around or bounce (but see 96] for problems 
with such bounces). 

Apart from negativity, the key difference between this "dark stiff matter" and the dark radiation term m/a 4 , is 
that the latter arises from the bulk Weyl curvature via the £ n V tensor, while the former arises from non- vacuum 



stresses in the bulk via the T^v tensor in Eq. ()3.2(J)) . 
modes of the graviton. 



The dark stiff matter does not arise from massive KK 
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Another example is provided by the Vaidya-AdSs metric, which can be written after transforming to a new 
coordinate v — T + J dR/F in Eq. (|5.1|) . so that v = const are null surfaces, and 
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(5.49) 
(5.50) 



This model has a moving FRW brane in a 4-isotropic bulk (which is not static), with either a radiating bulk 
black hole (dm/dv < 0), or a radiating brane (dm/dv > 0) [93. The metric satisfies the 5D field equations 1)3. 1|) 
with a null-radiation energy-momentum tensor, 



(5y T AB = tpk A k B , k A k A = , k A u J 



where tp cx dm/dv. It follows that 



(5.51) 



(5.52) 



In this case, the same effect, i.e. a varying mass parameter to, contributes to both £^ u and F^ v in the brane 
field equations. The modified Fricdmann equation has the standard 1-brane RS-type form, but with a dark 
radiation term that no longer behaves strictly like radiation: 
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By Eqs. (|3.25() and l|5.51|l . we arrive at the matter conservation equations, 



(5.53) 



(5.54) 



This shows how the brane loses (ip > 0) or gains (ip < 0) energy in exchange with the bulk black hole. For an 
FRW brane, this equation reduces to 



p + 3H(p + p) = -2ip. 
The evolution of m is governed by the 4D contracted Bianchi identity, using Eq. 1)5.52(1 : 



For an FRW brane, this yields 



where ps = 3m(i)/(re 2 a ). 
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(5.55) 
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• A more complicated bulk metric arises when there is a self-interacting scalar field $ in the bulk |29j, |40j • In the 
simplest case, when there is no coupling between the bulk field and brane matter, this gives 



{5) TAB = B ~ i5) 9AB 



(5.58) 



where $(x, y) satisfies the 5D Klein-Gordon equation, 

(5)D$_y'($) = 0. (5.59) 
The junction conditions on the field imply that 

d y ${x, 0) = . (5.60) 
Then Eqs. (|3.25|) and (|5.58|) show that matter conservation continues to hold on the brane in this simple case: 

VT pi/ = 0. (5.61) 
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From Eq. I|5.58[l one finds that 



where 



(x) = $(ac,0) 



so that the modified Friedmann equation becomes 
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When there is coupling between brane matter and the bulk scalar field, then the Friedmann and conservation 
equations are more complicated [2^, El . 



VI. BRANE- WORLD COSMOLOGY: PERTURBATIONS 



The background dynamics of brane-world cosmology are simple because the FRW symmetries simplify the bulk 
and rule out nonlocal effects. But perturbations on the brane immediately release the nonlocal KK modes. Then the 
5D bulk perturbation equations must be solved in order to solve for perturbations on the brane. These 5D equations 
are partial differential equations for the 3D Fourier modes, with both initial and boundary conditions needed. 

The theory of gauge- in variant perturbations in bran e -wor l d cos mology has been extensively investigated and de- 
veloped [H El lHllzl [zi HI |H HI |il H EE EH Ell and is qualitatively well understood. The key 
remaining task is integration of the coupled bran e-bu l k pe rturbation equations. Special cases have been solved, 
where the s e equations effectively d ecouple |68l 16^ . llOll Il02j , and approximation schemes have recently been devel- 
oped |l03l Il04 Il05l Il06t llOTt llOSfl for the more general cases where the coupled system must be solved. From the 
brane viewpoint, the bulk effects, i.e., the high-energy corrections and the KK modes, act as source terms for the 
brane perturbation equations. At the same time, perturbations of matter on the brane can generate KK modes (i.e., 
emit 5D gravitons into the bulk) which propagate in the bulk and can subsequently interact with the brane. This 
nonlocal interaction amongst the perturbations is at the core of the complexity of the problem. It can be elegantly 
expressed via integro-differential equations j^] , which take the form (assuming no incoming 5D gravitational waves) 

A k (t)= f dt'g(t,t')B k (t'), (6.1) 

where Q is the bulk retarded Green's function evaluated on the brane, and A kl B k are made up of brane metric and 
matter perturbations and their (brane) derivatives, and include high-energy corrections to the background dynamics. 
Solving for the bulk Green's function, which then determines G, is the core of the 5D problem. 

We can isolate the KK anisotropic stress 7r^ as the term that must be determined from 5D equations. Once 
is determined in this way, the perturbation equations on the brane form a closed system. The solution will be of the 
form, expressed in Fourier modes: 

Trf (t) oc fdt' Q{t,t')F k {t') , (6.2) 

where the functional F k will be determined by the covariant brane perturbation quantities and their derivatives. It is 
known in the case of a Minkowski background [3(]], but not in the cosmological case. 

The KK terms act as source terms modifying the standard general relativity perturbation equations, together with 
the high-energy corrections. For example, the linearization of the shear propagation equation l|3.82|l yields 

— * p 

&n„ + 2Ha fllJ + E^ v - y7i> - V ^A v) = yTrjL ~ ~( l + 3^)^ "> • ( 6 - 3 ) 

In 4D general relativity, the right hand side is zero. In the brane-world, the first source term on the right is the 
KK term, the second term is the high-energy modification. The other modification is a straightforward high-energy 
correction of the background quantities H and p via the modified Friedmann equations. 

As in 4D general relativity, there are various different, but essentially equivalent, ways to formulate linear cosmo- 
logical perturbation theory. First I describe the covariant brane-based approach. 
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A. 1+3-covariant perturbation equations on the brane 

In the 1+3-covariant approach |35lll0lifl09| . perturbative quantities are projected vectors, — V^) , and projected 
symmetric tracefree tensors, W^u = W/^ v \ , which are gauge-invariant since they vanish in the background. These are 
decomposed into (3D) scalar, vector and tensor modes as 

Vp. = + (6.4) 

= V {ll V v) W + V {li W v) +W^ v , (6.5) 

where — Wi^s and an overbar denotes a (3D) transverse quantity, 

V*% = = V"TV . (6.6) 

In a local inertial frame comoving with i// 1 , i.e., u M = (1,0), all time components may be set to zero: = (0, Vi), 

Purely scalar perturbations are characterized by the fact that vectors and tensors are derived from scalar potentials, 
i.e., 

Vju = = = . (6.7) 

Scalar perturbative quantities are formed from the potentials via the (3D) Laplacian, e.g., V = V M V M V^ = V 2 V '. 
Purely vector perturbations are characterized by 

V» = % , WV = V (tl W v) , curl V„/ = -2/^ , (6.8) 
where a> M is the vorticity, and purely tensor by 

V M / = = V„ , W„ v = W„ v ■ (6.9) 

The KK energy density produces a scalar mode V M p£ (which is present even if pg = in the background). The KK 
momentum density carries scalar and vector modes, and the KK anisotropic stress carries scalar, vector and tensor 
modes: 

4 = V M /+^, (6.10) 
- V^^+V^+fr^. (6.11) 

Linearizing the conservation equations for a single adiabatic fluid, and the nonlocal conservation equations, we 
obtain 

p + e(p + p) = 0, (6.12) 

(? s W tl + {p+p)A^ = Q, (6.13) 

ps + ^&Ps + V^=0, (6.14) 
q% + AHqi + )-V »p £ + \p £ A^ + V*V* = -^±^V M p . (6.15) 
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Linearizing the remaining propagation and constraint equations leads to 
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(6.16) 
(6.17) 
(6.18) 

(6.19) 

(6.20) 
(6.21) 
(6.22) 
(6.23) 
(6.24) 

(6.25) 



Equations l|6.12|l . I|6.14|l and l|6.16|l do not provide gauge-invariant equations for perturbed quantities, but their spatial 
gradients do. 

These equations are the basis for a 1+3-covariant analysis of cosmolo gical perturbations from the brane observer's 
viewpoint, following the approach developed in 4D general relativity |110| . The equations contain scalar, vector 
and tensor modes, which can be separated out if desired. They are not a closed system of equations until is 
determined by a 5D analysis of the bulk perturbations. An extension of the 1+3-covariant perturbation formalism 
to 1+4 dimensions would require a decomposition of the 5D geometric quantities along a timelike extension u A into 
the bulk of the brane 4-velocity field v?, and this remains to be done. The 1+3-covariant perturbation formalism is 
incomplete until such a 5D extension is performed. The metric-based approach does not have this drawback. 



Metric-based perturbations 



An alternative appr oach to brane-world cosmological perturba tion s is an extension of the 4D metric-based gauge- 
invariant theory 1 1 1 1| . A review of this approach is given in |l00j| . In an arbitrary gauge, and for a flat FRW 
background, the perturbed metric has the form 
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where the background metric functions A,N are given by Eqs. I|5.4|l and (|5.5|l . The scalars ip,"R., C,a, /3, v represent 
scalar perturbations. The vectors Si,Fi and Xi are transverse, so that they represent 3D vector perturbations, and 
the tensor fa is transverse traceless, representing 3D tensor perturbations. 

In the Gaussian normal gauge, the brane coordinate-position remains fixed under perturbation, 



^ds 2 = 



9uJ ( x > V) + Sg^{x, y) dx^dx v + dy z , 



(6.27) 



where gjS) is the background metric, Eq. (|5.3|l . In this gauge, 



a = j3 = v = Xi = 0- 



(6.28) 
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In the 5D longitudinal gauge, 



-B + C = = -(3 + C . 



(6.29) 



In this gauge, and for an AdS5 background, the metric perturbation quantities can all be expressed in terms of a 
"master variable" f2 which obeys a wave equation 98]. In the case of scalar perturbations, we have for example, 
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with similar expressions for the other quantities. All of the metric perturbation quantities are determined once a 
solution is found for the wave equation 
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The junction conditions Eq. (|3.19|) relate the off-brane derivatives of metric perturbations to the matter perturba- 
tions: 
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where 



ST° = -Sp, 5T\ = a 2 qi 
For scalar perturbations in the Gaussian normal gauge, this gives 
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where Sir is the scalar potential for the matter anisotropic stress, 



Sn,,i = didjSn - -S i3 d k d k 5n . 



The perturbed KK energy-momentum tensor on the brane is given by 

<5£% = K 2 5p £ , S£° l = -K 2 a 2 qf , 
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(6.36) 
(6.37) 
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(6.39) 

(6.40) 
(6.41) 



The evolution of the bulk metric perturbations is determined by the perturbed 5D field equations in the vacuum 
bulk, 



0. 



(6.42) 



Then the matter perturbations on the brane enter via the perturbed junction conditions, Eq. H6.32|) . 
For example, for scalar perturbations in Gaussian normal gauge, 
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For tensor perturbations (in any gauge) , the only nonzero components of the perturbed Einstein tensor are 
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In the following, I will discuss various perturbation problems, using either a 1+3-covariant or a metric-based 
approach. 

C. Density perturbations on large scales 

In the covariant approach, we define matter density and expansion (velocity) perturbation scalars, as in 4D general 
relativity, 



A = — V 2 p, Z 
P 

Then we can define dimensionless KK perturbation scalars [35j| . 



a 2 V 2 6. 



(6.45) 
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where the scalar potentials q and it are defined by Eqs. (|6.10|l and (|6.11() . The KK energy density (dark radiation) 
produces a scalar fluctuation U which is present even if pg = in the background, and which leads to a non-adiabatic 
(or isocurvature) mode, even when the matter perturbations are assumed adiabatic • We define the total effective 
dimensionless entropy Stot via 



Ptot S tat = a 2 V 2 ptot 
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where c 2 ot = ptot / ptot is given in Eq. (|3.6Q(I . Then 

9 k 2 - i + a + 



f) P/A] 



[3(1 + «;)(! + p/A) + 4p £ /p][3w + 3(1 + 2w)p/2A + p £ /p] 



4p£ 
3 P 



{l + w)U 



(6.47) 



(6.48) 



If pg — in the background, then f7 is an isocurvature mode: Stot oc {\-\-w)U. This isocurvature mode is suppressed 
during slow-roll inflation, when 1 + w ~ 0. 

If pg ^ in the background, then the weighted difference between U and A determines the isocurvature mode: 
Stot x (4pf /3p)A — (1 + w)U. At very high energies, p> A, the entropy is suppressed by the factor X/p. 

The density perturbation equations on the brane are derived by taking the spatial gradients of Eqs. I|6.12|l . (|6.14|) 
and lEIinj), and using Eqs. l|?TT3)) and !pTT5|) . This leads to H2 
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(6.49) 
(6.50) 

(6.51) 

(6.52) 



The KK anisotropic stress term II occurs only via its Laplacian, V 2 II. If we can neglect this term on large scales, 
then the system of density perturbation equations closes on super-Hubble scales |35j . An equivalent statement applies 
to the large-scale curvature perturbations [6^ . KK effects then introduce two new isocurv ature modes on large scales 
(associated with U and Q) , as well as modifying the evolution of the adiabatic modes [83, llOlj • 
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FIG. 8: The evolution of the covariant variable 3>, defined in Eq. I|6.55^ (and not to be confused with the Bardeen potential), 
along a fundamental world-line. This is a mode that is well beyond the Hubble horizon at N = 0, about 50 e-folds before 
inflation ends, and remains super-Hubble through the radiation era. A smooth transition from inflation to radiation is modelled 
by w = i[(2 — |e) tanh(iV — 50) — (1 — § e )], where e is a small positive parameter (chosen as e = 0.1 in the plot). Labels on 
the curves indicate the value of po/A, so that the general relativistic solution is the dashed curve (po/A = 0). (From 82].) 



Thus on large scales the system of brane equations is closed, and we can determine the density perturbations 
without solving for the bulk metric perturbations. 

We can simplify the system as follows. The 3-Ricci tensor defined in Eq. I|3.91|l leads to a scalar covariant curvature 
perturbation variable, 



C = a 4 \/ 2 R^ = -Aa 2 HZ + 2n 2 a 2 p (l + -?-) A + 2n 2 a 2 pU . 

V 2 A / 

It follows that C is locally conserved (along u M flow lines): 

C = C'o , Co = . 
We can further simplify the system of equations via the variable 

$ = n 2 a 2 pA . 



(6.53) 



(6.54) 



(6.55) 



This should not be confused with the Bardeen metric perturbation variable although it is the covariant analogue 
of $h in the general relativity li mit. I n the brane-world, high-energy and KK effects mean that <&h is a complicated 
generalization of this expression |l0l| , involving II, but the simple $ above is still useful to simplify the system of 
equations. Using these new variables, we find the closed system for large-scale perturbations: 



$ = —H 
U = -H 



l + (l + w) 
1 - 3w 
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2k 2 P£ 



(1 + w) 



2 4 2 

a K p 



3H 2 



U 



2p £ 
3a 2 Hp 



2H 
6c 2 H 2 



U 



{l + w) w 



A (1 + w)k 2 p 



Ps 



3a 2 Hp 



Co , 
Co- 



If there is no dark radiation in the background, ps = 0, then 

U = U Q exp - J {I - 3w)dN , 



(6.56) 
(6.57) 

(6.58) 
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and the above system reduces to a single equation for At low energies, and for constant w, the non-decaying 
attractor is the general relativity solution, 



At very high energies, for w > — i, 



_ 3(1 + w) 
$low ^2(5 + 3 U ,) C °' 



(6.59) 



$hi g h-> -—(1 + w) 
1 Po 



Co 2U 



7 + 6w 5 + 6u> 



(6.60) 



where U = h 2 o^pqUq, so that the isocurvature mode has an influence on Initially, <& is suppressed by the factor 
X/ po, but then it grows, eventually reaching the attractor value in Eq. (|6.59|) . For slow-roll inflation, when 1 + w ~ e, 
with < e < 1 and H~ l \e\ = \e'\ < 1, 



^high ~ -ze — Ooe 
* Po 



(6.61) 



where N = ln(a/ao), so that $ has a growing- mode in the early universe. This is different from general relativity, 
where $ is constant during slow-roll inflation. Thus more amplification of $ can be achieved than in general relativity, 
as discussed above. This is illustrated for a toy model of inflation-to-radiation in Fig. 8. The early (growing) and late 
time (constant) attractor solutions are seen explicitly in the plots. 
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FIG. 9: The evolution of <3> in the radiation era, with dark radiation present in the background. (From jl IS 



The presence of dark radiation in the b ackgr ound introduces new features. In the radiation era (w — h), the 
non-decaying low-energy attractor becomes |l!2| 

$iow « ^y(l-a), (6.62) 

a = — ^ 0.05 . (6.63) 



41 



The dark radiation serves to reduce the final value of $, leaving an imprint on <I>, unlike the ps 
In the very high energy limit, 



case, Eq. <|H35|l . 



high 



po 



Co 



16a 



273 



4^o 
539 



(6.64) 



Thus $ is initially suppressed, then begins to grow, as in the no-dark-radiation case, eventually reaching an attractor 
which is less than the no-dark-radiation attractor. This is confirmed by the numerical integration shown in Fig. 9. 



D. Curvature perturbations and the Sachs- Wolfe effect 

The curvature perturbation 1Z on uniform density surfaces is defined in Eq. (|6.26|) . The associated gauge-invariant 
quantity 

< =n+ w^y (6 - 65) 

may be defined for matter on the brane. Similarly, for the Weyl "fluid" if pg ^ in the background, the curvature 
perturbation on hypersurfaces of uniform dark energy density is 

C £ = K + 5 -p- . (6.66) 

On large scales, the perturbed dark energy conservation equation is |6S| 

(S P£ y + W5p £ + 4:p £ K = 0, (6.67) 

which leads to 

6 = 0. (6.68) 
For adiabatic matter perturbations, by the perturbed matter energy conservation equation, 

(dp)' + 3H{Sp + Sp) + 3(p + p)U = , (6.69) 

we find 

C = 0. (6.70) 

This is independent of brane- world modifications to the field equations, since it depends on energy conservation only. 
For the total, effective fluid, the curvature perturbation is defined as follows [68| : if ps ^ in the background, 



Ctot — c + 

and if p£ = in the background, 



4p £ 



3(p + p)(l + p/\)+4p £ 



(& - , (6.71) 



C " = C+ 3(P + ,KWA) (6 ' 72) 
S P£ = 6 -^f, (6.73) 

where SCs is a constant. It follows that the curvature perturbations on large scales, like the density perturbations, 
can be found on the brane without solving for the bulk metric perturbations. 

Note that ^tot 7^ even for adiabatic matter perturbations; for example, if ps = in the background, 

The KK effects on the brane contribute a non-adiabatic mode, although Ctot — > at low energies. 
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Although the density and curvature perturbations can be found on super-Hubble scales, the Sachs- Wolfe effect 
requires in order to translate from density/ curvature to metric perturbations. In the 4D longitudinal gauge of 
the metric perturbation formalism, the gauge-invariant curvature and metric perturbations on large scales are related 
by 

(6.75) 
(6.76) 

where the radiation anisotropic stress on large scales is neglected, as in general relativity, and Sir £ is the scalar 
potential fo r tt^ v , equivalent to the covariant quantity II defined in Eq. I|6.46[) . In 4D general relativity, the right hand 
side of Eq. (|6.76[) is zero. The (non-integrated) Sachs- Wolfe formula has the same form as in general relativity: 




ST 
T 



= (Crad+V>-ft)|dec- (6.77) 

The brane-world corrections to the general relativistic Sachs- Wolfe effect are then given by j6^| 

ST ( ST\ 8 / p rad 



T \ T J 3 \Pcdr 



S £ - K 2 a 2 Sn £ + / da a 7/2 Stt £ , (6.78) 
opl z J 



where S £ is the KK entropy perturbation (determined by 8p £ ). The KK term &ix £ cannot be determined by the 4D 
brane equations, so that ST /T cannot be evaluated on large scales without solving the 5D equations. [Equation .7811 
has been generalized to a 2-brane model, in which the radion makes a contribution to the Sachs- Wolfe effect |l!3j |.] 
The presence of the KK (Weyl, dark) component has essentially two possible effects. 

• A contribution from the KK entropy perturbation Sg that is similar to an extra isocurvature contribution. 

• The KK anisotropic stress Sir £ also contributes to the CMB anisotropies. In the absence of anisotropic stresses, 
the curvature perturbation (tot would be sufficient to determine the metric perturbation 1Z and hence the large- 
angle CMB anisotropies, via Eqs. (|6.75|) . I|6.76|l and l|6.77|l . However bulk gravitons generate anisotropic stresses 
which, although they do not affect the large-scale curvature perturbation (tot, can affect the relation between 
(tot, TZ and ip, and hence can affect the CMB anisotropies at large angles. 

A simple phenomenological approximation to 5ir £ on large scales is discussed in j6^| and the Sachs- Wolfe effect is 
estimated as 



T V 9 /in \idec 



1 '■ r ln(/ h /; 



ln(i oq /^ 



(6.79) 



where £4 is the 4D Planck time and t- ln is the time when the KK anisotropic stress is induced on the brane, which is 
expected to be of the order of the 5D Planck time. 

A self-consistent approximation is developed in |104| . using the low-energy 2-brane approximation |103| to find an 
effective 4D form for and hence for 8tt £ . This is discussed below. 

E. Vector perturbations 

The vorticity propagation equation on the brane is the same as in general relativity, 



u tl +2Hu tl = --cm\A, x . (6.80) 



Taking the curl of the conservation equation l|3.63|) (for the case of a perfect fluid, = = n^), and using the 
identity in Eq. Q6.8[l . one obtains 

curl^ = -6£Tc^ , (6.81) 
as in general relativity, so that Eq. (|6.80(l becomes 

ci„ + (2 - 3c 2 ) ffw M = , (6.82) 
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which expresses the conservation of angular momentum. In general relativity, vector perturbations vanish when the 
vorticity is zero. By contrast, in brane- world cosmology, bulk KK effects can source vector perturbations even in the 
absence of vorticity ^4|. This can be seen via the divergence equation for the magnetic part of the 4D Weyl 
tensor on the brane: 



V 2 H^ = 2k 2 ( P + P ) 



-K psuj^ - -k curl <7 



3 2' 



(6.83) 



where H^ v = S/^H^y Even when to^ = 0, there is a source for gravimagnetic terms on the brane from the KK 
quantity curl . 

We define covariant dimensionless vector perturbation quantities for the vorticity and the KK gravi-vector term: 



= ftp = -curlgj 



On large scales, we can find a closed system for these vector perturbations on the brane [14j: 

a M + (1 - 3c s 2 ) JTa„ = 0, 



4(3^-1)^-9(1 



w) - 



(6.84) 

(6.85) 
(6.86) 



Thus we can solve for and /3 M on super-Hubble scales, as for density perturbations. Vorticity in the brane matter 
is a source for the KK vector perturbation /3 M on large scales. Vorticity decays unless the matter is ultra-relativistic or 
stiffcr (w > and this source term typically provides a decaying mode. There is another pure KK mode, independent 
of vorticity, but this mode decays like vorticity. For w = pj p = const, the solutions are 





f a 












f a 











3w-l 



opo \a J A \ a 



(6.87) 
(6.88) 



where = = c M . 

Inflation will redshift a way the vorticity and the KK mode. Indeed, the massive KK vector modes are not excited 
during slow-roll inflation |l00j ] . 



F. Tensor perturbations 

The covariant description of tensor modes on the brane is via the shear, which satisfies the wave equation ^ 



2A+i K 2 {p-S P -(p + 3p)j} 



(6.89) 



Unlike the density and vector perturbations, there is no closed system on the brane for large scales. The KK anisotropic 
stress is an unavoidable source for tensor modes on the brane. Thus it is necessary to use the 5D metric-based 
formalism. This is the subject of the next section. 

VII. GRAVITATIONAL WAVE PERTURBATIONS IN BRANE- WORLD COSMOLOGY 



The tensor perturbations are given by Eq. I|6.26[l . i.e. (for a flat background brane), 

< 5 W = -N 2 (t, y)dt 2 + A 2 {t, y) [<% + dx l dx j + dy 2 . 



(7.1) 



The transverse traceless satisfies Eq. 

f(t,v), 



1 

V 2 " 



which implies, on splitting into Fourier modes with amplitude 
A' N n 



A* f = f " + [ 3 'A + N^' 



(7.2) 
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By the transverse traceless part of Eq. (|6.32(l , the boundary condition is 



f 



ij Ibrano 



(7.3) 



where 7fjj is the tensor part of the anisotropic stress of matter-radiation on the brane. 

The wave equation 1|7.2[1 cannot be solved analytically except if the background metric functions are separable, and 
this only happens for maximally symmetric branes, i.e. branes with constant Hubble rate H . This includes the RS 
case H = 0, already treated in Sec. II. The cosmologically relevant case is the de Sitter bran e, H > 0. We can 
calculate the spectrum of gravitational waves generated during brane inflation [g^, EE Eill If , if we approximate 
slow-roll inflation by a succession of de Sitter phases. The metric for a de Sitter brane, dS4, in AdSs is given by 
Eqs. (E31)-(E21) with 



N(t,y) 

n(y) = cosh/iy — ( 1 



n (y), A (t, V) = a{t)n{y) . 
A 



sinh^|y| 



i(t) = do exp Ho(t — to) . 



M 2 - K n (\ -X- P ° 



(7.4) 
(7.5) 

(7.6) 



where /x 





FIG. 10: Graviton "volcano" potential around the dS4 brane, showing the mass gap (from [Tfl p 



The linearized wave equation l|7.2|l is separable. As before, we separate the amplitude as / = ^ <Pm{t) f m {y) where 
m is the 4D mass, and this leads to: 



ifim + 3H <p m + 



k 2 



a* 



ip m = 0. 



.C + 4-C + ^/ m = o. 

n n 



The general solutions for m > are 



r -,„(/ ! = exp j --//,,/ ) U,, ( — e 



-Hot 



f m (y) = n {y)-^Ll /2 [ x jl+^n{y) 



(7.7) 
(7.8) 

(7.9) 
(7.10) 



where B v is a linear combination of Bessel functions, L^^ 2 is a linear combination of associated Legendre functions, 
and 
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It is more useful to reformulate Eq. (|7.8(l as a Schrodinger-type equation, 

d 2 <Z> r 



dz 2 



V(z)* r , 



using the conformal coordinate 



z = z h 



v dy 1 . ! /JJo 

-T^r > 2 b = — smh 

o n(y) H \ \x 



and defining ^ m = n 3 / 2 f m . The potential is given by (see Fig. 10) 

9 



V(z) 



15H 2 



4sinh 2 (Jf z) 4 



H 2 - 3/i 



1 



Po 
A 



6(z - Zb) , 



(7.12) 



(7.13) 



(7.14) 



where the last term incorporates the boundary condition at the brane. The "volcano" shape of the potential show s 
how the 5D graviton is localized at the brane at low energies. (Note that localization fails for an AdS4 brane |115| .') 
The non-zero value of the Hubble parameter implies the existence of a mass gap 26], 



Am = - Hq , 



(7.15) 



between the zero mode and the continuum of massive KK modes. This result has been generalized: for dS4 brane(s) 
with bulk scalar field, a universal lower bound on the mass gap of the KK tower is js3 



(7.16) 



Am > \l-Ho. 



The massive modes decay during inflation, according to Eq. (|7.9|l . leaving only the zero mode, which is effectively a 
4D gravitational wave. The zero mode, satisfying the boundary condition, 



is given by 



where the normalization condition 



implies that the function F is given by [8f 



/o(^0)=0, 
fa = sfji F(H /n) , 

2 |¥g|d*=l, 



(7.17) 



(7.18) 



(7.19) 



F(x) = I \j\ + x 2 - x 2 In 



,-1/2 



(7.20) 



At low energies {Ho -C fj), we recover the general relativity amplitude: F — * 1. At high energies, the amplitude is 
considerably enhanced: 



H > fi =*> F 



I 3H Q 
2/x 



(7.21) 



The factor F determines the modification of the gravitational wave amplitude relative to the standard 4D result: 

2" 



A 



Hi 

M 2 \ 2tt 



F 2 (H /fi). 



(7.22) 



The modifying factor F can also be interpreted as change in the effective Planck mass |83j. 
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FIG. 11: Damping of brane-world gravity waves on horizon re-entry due to massive mode generation. The solid curve is the 
numerical solution, the short-dashed curve the low-energy approximation, and the long-dashed cur ve t he standard general 
relativity solution, eg = po/X and 7 is a parameter giving the location of the regulator brane. (From |106(I .') 



This enhanced zero mode produced by brane inflation remains frozen outside the Hubble radius, as in general 
relativity, but when it re-enters the Hubble radius during radiation or matter domination, it will no longer be separated 
from the massive modes, since H will not be constant. Instead, massive modes will be excited during re-entry. In 
other words, energy will be lost from the zero mode as 5D gravitons are emitted into the bulk, i.e., as massive modes 
are produc ed on the brane. A phenomenological model of the damping of the zero mode due to 5D graviton e mission 
is given in |l02j | . Self-consistent low-energy approximations to compute th is eff ect are developed in |l06t Il07| . 

At zero order, the low-energy approximation is based on the following |108| . In the radiation era, at low energy, 
the background metric functions obey 

A(t, y) -> a(t)e-"» , N(t, y) -» e"™ . (7.23) 

To lowest order, the wave equation therefore separates, and the mode functions can be found analytically 108]. The 
massive modes in the bulk, f m (y), are the same as for a Minkowski brane. On large scales, or at late times, the mode 
functions on the brane are given in conformal time by 

^(V) = V- 1/2 B 1/ J^V 2 ) , (7.24) 



V2 



where ah marks the start of the low-energy regime (ph = A) and B v denotes a linear combination of Bessel functions. 
The massive modes decay on super- Hubble scales, unlike the zero- mode. Expanding the wave equation in po/X, one 

arrives at the first order, where mode-mixing arises. The massiv e modes ipl^ (ry) on sub-Hubble scales are sourced by 
the initial zero mode that is re-entering the Hubble radius 107]: 



atp$ + k 2 a<p£> + m 2 a 3 (p$ = -4^/ m0 fc 2 a^ 0) , (7.25) 



where I m o is a transfer matrix coefficient. The numerical integration of the equations 106] confirms the effect of 
massive mode generation and consequent damping of the zero- mode, as shown in Fig. 11. 
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VIII. CMB ANISOTROPIES IN BRANE- WORLD COSMOLOGY 



For the CMB anisotropies, one needs to consider a multi-component source, 
expressions for the total effective energy-momentum tensor, we obtain 



Linearizing the general nonlinear 



Ptot 

Ptot 
tot 

~tot 



Mi 



p 



2A 



P_ 

2A 

(2p- 



1 



P 

P)- 
h3 P 



P£ 

3 



2A 



(8.1) 
(8.2) 
(8.3) 
(8.4) 



where 



) . 9m 



(0 
A 1 



■5) 



are the total matter-radiation density, pressure and momentum density, and 71"^,, is the photon anisotropic stress 
(neglecting that of neutrinos, baryons and CDM). 

The perturbation equations in the previous section form the basis for an analysis of scalar and tensor CMB 
anisotropies in the brane-world. The f ull system o f equ ations on the brane, including the Boltzmann equation for 
photons, has been given for scalar |lOl| and tensor |l02fl perturbations. But the systems are not closed, as discussed 
above, because of the presence of the KK anisotropic stress 7r^ v , which acts a source term. 

In the tight-coupling radiation era, the scalar perturbation equations may be decoupled to give an equation for the 
gravitational potential defined by the electric part of the brane Weyl tensor (not to be confused with £ M „): 



(8.6) 



In general relativity, the equation in <E> has no source term, but in the brane-world there is a source term made up of 
tx^ v and its time-derivatives. At low energies (p <C A), and for a flat background (K = 0), the equation is |l0l| 
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- 7T 
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where x = k/(aH), a prime denotes d/dx, and $>k, Trf are the Fourier modes of $ and tt^. In general relativity the 
right hand side is zero, so that the equation may be solved for <I>/j, and then for the remaining perturbative variables, 
which gives the basis for initializing CMB numerical integrations. At high energies, earlier in the radiation era, the 
decoupled equation is fourth order 101]: 
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(8.8) 



The formalism and machinery are ready to compute the temperature and polarization anisotropies in brane-world 
cosmology, once a solution, or at least an approximation, is given for 7i~„. The resulting power spectra will reveal 
the nature of the brane-world imprint on CMB anisotropics, and would in principle provide a means of constraining 
or possibly falsifying the brane-world models. Once this is achieved, the implications for the fundamental underlying 
theory, i.e. M theory, would need to be explored. 

However, the first step required is the solution for ~k„ v . This solution will be of the form given in Eq. |JB2Jl. Once 
Q and F k are determined or estimated, the numerical integration in Eq. I|6.2|) can in principle be incorporated into a 
modified version of a CMB numerical code. The full solution in this form represents a formidable problem, and one 
is led to look for approximations. 
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A. The low-energy approximation 



The basic idea of the low-energy approximation 103J is to use a gradient expansion to exploit the fact that, during 
most of the history of the universe, the curvature scale on the observable brane is much greater than the curvature 
scale of the bulk [I < 1 mm): 



l ~ |iWT 1/2 » e~ I (5) Rabcd\ 
=> |v M | ~i- x «: \a y \ 



-1/2 



These conditions are equivalent to the low energy regime, since £ 2 oc A 1 and |i?^ va; a| ~ I^m^I' 
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L 2 



A 



< 1 



(8.9) 



(8.10) 



Using Eq. I|8.9|l to neglect appropriate gradient terms in an expansion in £ 2 /L 2 , the low-energy equations can be 
solved. However, two boundary conditions are needed to determine all functions of integration. This is achieved by 
introducing a second brane, as in the RS 2-brane scenario. This brane is to be thought of either as a regulator brane, 
whose backreaction on the observable brane is neglected (which will only be true for a limited time), or as a shadow 
brane with physical fields, which have a gravitational effect on the observable brane. 

The background is given by low-energy FRW branes with tensions ±A, proper times t±, scale factors a±, energy 
densities p± and pressures p± , and dark radiation densities ps ± . The physical distance between the branes is £d{t) , 
and 
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dt 



— = e a — , a_ = a+e- d , H_=e a [H+-d 
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(8.11) 



Then the background dynamics is given by (see [ill 

for the general background, including the high-energy regime): 

(8.12) 



H 2 ± =±—(p ± ± P£± ) 
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.13) 



The dark energy obeys pg+ = C/a 4 ^, where C is a constant. From now on, we drop the -{—subscripts which refer to 
the physical, observed quantities. 

The perturbed metric on the observable (positive tension) brane is described, in longitudinal gauge, by the metric 
perturbations tp and 1Z, and the perturbed radion is d = d + N. The approximation for the KK (Weyl) energy- 
momentum tensor on the observable brane is 



ft* 



„2d 



-y(^ + e- 2d T%) 



1 



and the field equations on the observable brane can be written in scalar-tensor form as 

G\ = —TZ+ K2{1 ~ X)2 T» V 

X X 



.14) 



1 



.15) 



where 



X = 1 



-2d 



21-X 



(8.16) 



The perturbation equations can then be derived as generalizations of the standard equations. For example, the 
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8G° equation is |l!7| 
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.17) 



The trace part of the perturbed field equation shows that the radion perturbation determines the crucial quantity, 
Sir £ : 



K + ip = — N = -K 2 a 2 Sir £ , 

e 2d _ I 

where the last equality follows from Eq. (|6.76(l . The radion perturbation itself satisfies the wave equation 



(8.18) 
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A new set of variables (p± , E turns out be very useful |l04l Ill3j | : 
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.20) 



.21) 



The variable E determines the metric shear in the bulk, whereas tp± give the brane displacements, in transverse 
traceless gauge. The latter variables have a simple relation to the curvature perturbations on large scales |104 ll 13) 
(restoring the +-subscripts): 



Ctot ± 



-tf± 



.22) 



where f± = df±/dt±. 



B. The simplest model 

The simplest model is the one in which 

p £ = = d (8.23) 

in the background, with p-/ p- = p/p. The regulator brane is assumed to be far enough away that its effects on the 
physical brane can be neglected over the timescales of interest. By Eq. Q8.13p . it follows that 

P-=-pe 2d , (8.24) 
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i.e., the matter on the regulator brane must have fine-tuned and negative energy density to prevent the regulator brane 
from moving in the background. With these assumptions, and further assuming adiabatic perturbations for the matter, 
there is only one independent brane- world parameter, i.e., the parameter measuring dark radiation fluctuations: 
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FIG. 12: The CMB power spectrum with brane-world effects, encoded in the dark radiation flu ctu ation parameter <5C» as a 
proportion of the large-scale curvature perturbation for matter (denoted in the plot). (From [1Q4|.) 



This assumption has a remarkable consequence on large scales: the Weyl anisotropic stress Sirg terms in the Sachs- 
Wolfe formula Eq. I|tj.78fl cancel the entropy perturbation from dark radiation fluctuations, so that there is no difference 
on the largest scales from the standard general relativity power spectrum. On small scales, beyond the first acoustic 
peak, the brane-world corrections are negligible. On scales up to the first acoustic peak, brane-world effects can be 
significant, changing the height and the location of the first peak. These features are apparent in Fig. 12. However, it 
is not clear to what extent these features are general brane-world features (within the low-energy approximation), and 
to what extent they are consequences of the simple assumptions imposed on the background. Further work remains 
to be done. 

A related low-energy approximation, using the moduli space approximation, has been developed for certain 2-brane 
models with bulk scalar field |105|| . The effective gravitational action on the physical brane, in the Einstein frame, is 



Sefi = 



1 



d 4 Xy 



R 



12a 2 
l + 2a 2 



W) 2 



G 



l + 2a 2 



(8.26) 



where a is a coupling constant, and <j> and x are moduli fields (determined by the zero-mode of the bulk scalar field 
and the radion). Figure 13 shows how the CMB anisotropies are affected by the x-field. 



IX. CONCLUSION 



Simple brane-world models of RS type provide a rich phenomenology for exploring some of the ideas that are 
emerging from M theory. The higher-dimensional degrees of freedom for the gravitational field, and the confinement 
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of standard model fields to the visible brane, lead to a complex but fascinating interplay between gravity, particle 
physics and geometry, that enlarges and enriches general relativity in the direction of a quantum gravity theory. 

This review has attempted to show some of the key features of brane-world gravity from the perspective of astro- 
physics and cosmology, emphasizing a geometric approach to dynamics and perturbations. It has focused on 1-brane 
RS-typc brane- worlds which have some attractive features: 

• they provide a simple 5D phenomenological realization of the Horava-Witten supergravity solutions in the limit 
where the hidden brane is removed to infinity and the moduli effects from the 6 further compact extra dimensions 
may be neglected; 

• they develop a new geometrical form of dimensional reduction based on a strongly curved (rather than flat) 
extra dimension; 

• they provide a realization to lowest order of the AdS/CFT correspondence; 

• they incorporate the self-gravity of the brane (via the brane tension); 

• they lead to cosmological models whose background dynamics are completely understood and reproduce general 
relativity results with suitable restrictions on parameters. 

The review has highlighted both the successes and the remaining open problems of the RS models and their 
generalizations. The open problems stem from a common basic difficulty, i.e., understanding and solving for the 
gravitational interaction between the bulk and the brane (which is nonlocal from the brane viewpoint). The key open 
problems of relevance to astrophysics and cosmology are: 

• to find the simplest realistic solution (or approximation to it) for an astrophysical black hole on the brane, and 
settle the questions about its staticity, Hawking radiation and horizon; 

• to develop realistic approxim ation schemes (buil ding on recent work |l03l llo4. Il05l llOri Il07jp and manageable 
numerical codes (building on |!04ll05lll0fllTcn| ) to solve for the cosmological perturbations on all scales, com- 
pute the CMB anisotropics and large-scale structure, and impose observational constraints from high-precision 
data. 
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The RS-type models are the simplest brane- worlds with curved extra dimension that allow for a meaningful approach 
to astrophysics and cosmology. One also needs to consider generalizations that attempt to make these models more 
realistic, or that explore other aspects of higher-dimensional gravity which are not probed by these simple models. 
Two important types of generalization are the following. 

• The inclusion of dynamical interaction between the brane(s) and a bulk scalar field, so that the 
action is 



1 



dW- ( % \ {5) R ~ n 2 5 d A <Z>d A <5> - 2A 5 ($) 



d x-\/—g 



branc(s) 



-A($) 



K 



(9.1) 



(See jHIIlEt) 

The scalar field could represent a bulk dilaton of the gravitational sector, or a modulus field encoding the 
dynamica l influence on the effective 5D theory of an extra dimension other than the large fifth dimension |75| 
ESllii3. For two-b rane models, the brane separation introduces a new scalar degree of freedom, the radion. 
For general potentials of the scalar field which provide radion stabilization, 4D Einstein gravity is recovered 
at low energies on either brane p2^ . (By contrast, in the absence of a bulk scalar, low energy gravity is of 
Brans-Dicke type pfj.) 

In particular, such models will allow some fundamental problems to be addressed: 

— the hierarchy problem of particle physics; 

— an extra-dimensional mechanism for initiating inflation (or the hot radiation era with super-Hubble corre- 
lations) via brane interaction (building on the initial work in [l2> Iz3> 13 Iz3 ) i 

— an extra-dimensional explanation for the dark energy (and possibly also dark matter) puzzles: could dark 
energy or late-time acceleration of the universe be a result of gravitational effects on the visible brane of 
the shadow brane, mediated by the bulk scalar field? 

• The addition of stringy and quantum corrections to the Einstein-Hilbert action, including: 

— Higher-order curvature invariants, which arise in the AdS/CFT correspondence as next-to-leading order 
corrections in the CFT. The Gauss-Bonnet combination in particular has unique properties in 5D, giving 
field equations which are second-order in the bulk metric (and linear in the second derivatives), and being 
ghost-free. The action is 



S 



1 

24 
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d iv'-i 



-A + 



K 



(9.2) 



where a is the Gauss-Bonnet couplin g con stant, related to the string scale. The cosmological dynamics of 
these brane- worlds is investigated in |l!9| . In [5(j it is shown that the black string solution of the form of 
Eq. (|4.2|l is ruled out by the Gauss-Bonnet term. In this sense, the Gauss-Bonnet correction removes an 
unstable and singular solution. 

In the early universe, the Gauss-Bonnet corrections to the Friedmann equation have the dominant form 



H 2 oc p 2 / 3 : 



(9.3) 



at the highest energies. If the Gauss-Bonnet term is a small correction to the Einstein-Hilbert term, as 
may be expected if it is the first of a series of higher-order corrections, then there will be a regime of 
RS-dominance as the energy drops, when H 2 oc p 2 . Finally at energies well below the brane tension, the 
general relativity behaviour is recovered. 

Quantum field theory corrections arising from the coupling between brane matter and bulk gravitons, 
leading to an induced 4D Ricci term in the brane action. The original induced gravity brane-world |l20j 
was put forward as an alternative to the RS mechanism: the bulk is flat Minkowski 5D spacetime (and 
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as a consequence there is no normalizable zero- mode of the bulk graviton), and there is no brane tension. 
Another viewpoint is to see the induced-gravity term in the action as a correction to the RS action: 



where (3 is a positive coupling constant. Unlike the other brane- worlds discussed, these models lead to 5D 
behaviour on large scales rather than small scales. The cosmologi cal m odels have been analyzed in |l2l| . 
(Brane- world black holes with induced gravity are investigated in |122| .') 

The late-universe 5D behaviour of gravity can naturally produce a late-time acceleration, even without dark 
energy, although the fine-tuning problem is not evaded. 

The effect of the induced-gravity correction at early times is to restore the standard behaviour H 2 oc p to 
lowest order at the highest energies. As the energy drops, but is still above the brane tension, there is an 
RS regime, H 2 oc p 2 . In the late universe at low energies, instead of recovering general relativity, there are 
strong deviations from general relativity, and late-time acceleration from 5D gravity effects (rather than 
negative pressure energy) is typical. 

Thus we have a striking result that both forms of correction to the gravitational action, i.e., Gauss-Bonnet 
and induced gravity, suppress the Randall-Sundrum type high-energy modifications to the Friedmann 
equation when the energy reaches a critical level. ( Cosm ologies with both induced-gravity and Gauss- 
Bonnet corrections to the RS action are considered in |l23f .1 



In summary, brane-world gravity opens up exciting prospects for subjecting M theory ideas to the increasingly 
stringent tests provided by high-precision astronomical observations. At the same time, brane-world models provide 
a rich arena for probing the geometry and dynamics of the gravitational field and its interaction with matter. 
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